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Abstract 

We consider one-dimensional, locally finite interacting particle systems with two con- 
servation laws which under Eulerian hydrodynamic limit lead to two-by-two systems of 
conservation laws: 

f d tP + d x v{p,u) = o 

\ d t u + d x $(p,u) = 0, 

with (p, u) € T> C R 2 , where I? is a convex compact polygon in R 2 . The system is typically 
strictly hyperbolic in the interior of T> with possible non-hyperbolic degeneracies on the 
boundary dT>. We consider the case of isolated singular (i.e. non hyperbolic) point on the 
interior of one of the edges of T>, call it (po,Uo) = (0,0) and assume T> C {p > 0}. This 
can be achieved by a linear transformation of the conserved quantities. We investigate the 
propagation of small nonequilibrium perturbations of the steady state of the microscopic 
interacting particle system, corresponding to the densities (po, uq) of the conserved quanti- 
ties. We prove that for a very rich class of systems, under proper hydrodynamic limit the 
propagation of these small perturbations are universally driven by the two-by-two system 

f d t p + d x (pu) = 

\ d t u + d x {p + ~fu 2 ) = 

where the parameter 7 := ^^> uu (po, uo) (with a proper choice of space and time scale) is the 
only trace of the microscopic structure. The proof is valid for the cases with 7 > 1. 

The proof essentially relies on the relative entropy method and thus, it is valid only 
in the regime of smooth solutions of the pde. But there are essentially new elements: in 
order to control the fluctuations of the terms with Poissonian (rather than Gaussian) decay 
coming from the low density approximations we have to apply refined pde estimates. In 
particular Lax entropies of these pde systems play a not merely technical key role in the 
main part of the proof. 
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1 Introduction 

1.1 The PDE to be derived and some facts about it 

We consider the pde 

J d t p + d x (pu) = 

\ d t u + d x ( P + ju 2 ) =o 



i.i) 



where p = p(t,x) E [0, oo), u = u(t,x) S (—00,00) are density, respectively, velocity field and 
7 G R is a fixed parameter. For any fixed 7 this is a hyperbolic system of conservation laws in 
the domain (p, u) G M. + x R. 

Phenomenologically, the pde describes a deposition/ domain growth - or, in biological term: 
chemotaxis - mechanism: p(t, x) is the density of population performing the deposition and 
h(t, x) is the height of the deposition. Let 

u(t,x) := —d x h(t,x). 

The physics of the phenomenon is contained in the following two rules: 

(a) The velocity field of the population is proportional to the negative gradient of the height 
of the deposition. That is, the population is pushed towards the local decrease of the 
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deposition height. This rule, together with the conservation of total mass of the population 
leads to the continuity equation (the first equation in our system). 

(b) The deposition rate is 

d t h = p + -f(d x h) 2 . 

The first term on the right hand side is just saying that deposition is done additively by the 
population. The second term is a self-generating deposition, introduced and phenomeno- 
logically motivated by Kardar-Parisi-Zhang [5] and commonly accepted in the literature. 
Differentiating this last equation with respect to the space variable x results in the second 
equation of our system. 

The pde is invariant under the scaling: 

p(t, x) := A 2l3 p(A 1+ ^t, Ax), u(t, x) := A p u(A 1+f> t, Ax), 

where A > and (3 £ R are arbitrarily fixed. The choice (3 = gives the straightforward 
hyperbolic scale invariance, valid for any system of conservation laws. More interesting is the 
(3 = 1/2 case. This is the natural scale invariance of the system, since the physical variables 
(density and velocity fields) change covariantly under this scaling. This is the (presumed, 
but never rigorously proved) asymptotic scale invariance of the Kardar-Parisi-Zhang deposition 
phenomena. The nontrivial scale invariance of the pde Q1-1J) suggests its universality in some 
sense. Our main result indeed states its validity in a very wide context. 

It is also clear that the pde is invariant under the left-right reflection symmetry x i— ► — x: 
The parameter 7 of the pde (jl.lj) is of crucial importance: different values of 7 lead to 
completely different behavior. Here are listed some particular cases which arose in the past in 
various contexts: 

- The pde with 7 = arose in the context of the 'true self-repelling motion' constructed 
by Toth and Werner in j2B]- For a survey of this case see also [21]. The same equation, with 
viscosity terms added, appear in mathematical biology under the name of (negative) chemotaxis 
equations (see e.g. [Ej, [H]) ). 

- Taking 7 = 1/2 we get the 'shallow water equation'. See [3], US]- This is the only value of the 
parameter 7 when m = pu is conserved and as a consequence the pde 1)1- ljl can be interpreted 
as gas dynamics equation. 

- With 7 = 1 the pde is called 'Leroux's equation' which is of Temple class and for this 
reason much investigated. For many details about this equation see I n the recent paper 
[5] Leroux's system has been derived as hydrodynamic limit under Eulerian scaling for a two- 
component lattice gas, going even beyond the appearance of shocks. 

The main facts about the pde <fl.l|> are presented in Section ITT1 Here we only mention that 
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1. For any 7 G R the system is strictly hyperbolic in (p, u) £ (0, 00) x R, with hyperbolicity 
marginally lost at (p,u) = (0,0) for 7 / 1/2 and at p = for 7 = 1/2. This follows from 
straightforward computations. 

2. The Riemann invariants (or characteristic coordinates) are explicitly computed in section 
111! for a first impression see Figure |3] of the Appendix where the level lines of the Riemann 
invariants are shown. It turns out that the picture changes qualitatively at the critical values 
7 = 1/2, 7 = 3/4 and 7 = 1. It is of crucial importance for our later problem that the level 
curves, expressed as u 1— > p(u) are convex for 7 < 1, linear for 7 = 1 and concave for 7 > 1. 

3. For any 7 > the system (|l.l(l is genuinely nonlinear in (p, u) £ (0, 00) x R, with genuine 
nonlinearity marginally lost at (p, u) = (0, 0) for 7 7^ 0, 1/2 and at p = for 7 = 0, 1/2. (For 
7 < genuine nonlinearity is lost on the parabola p = —47(27 — 1) 2 (7 + l)~ 2 u 2 .) 

4. The system is sufficiently rich in Lax entropies. 

5. For 7 > the system I)1.1JI satisfies the conditions of the Lax-Chuey-Conley-Smoller Maxi- 
mum Principle (see (T^l, 

From the Maximum Principle a very essential difference between the cases 7 < 1,7 = 1 and 
7 > 1 follows, which is of crucial importance for our further work. In the case 7 < 1 all convex 
domains bounded by level curves of the Riemann invariants are unbounded (non-compact) and 
thus there is no a priori bound on the entropy solutions. Even starting with smooth initial data 
with compact support nothing prevents the solutions to blow up indefinitely. On the other hand, 
if 7 > 1 any bounded subset of R+ x R is contained in a compact convex domain bounded by 
level sets of the Riemann invariants, which fact yields a priori bounds on the entropy solutions, 
given bounded initial data. 

The goal of the present paper is to derive the two-by-two hyperbolic system of conservation 
laws (|1.1|) as decent hydrodynamic limit of some systems of interacting particles with two 
conserved quantities. 

We consider one-dimensional, locally finite interacting particle systems with two conserva- 
tion laws which under Eulerian hydrodynamic limit lead to two-by-two systems of conservation 
laws 



with (p,u) £ V C R 2 , where V is a convex compact polygon in R . The system is typically 
strictly hyperbolic in the interior of T> with possible non-hyperbolic degeneracies on the bound- 
ary dT>. We consider the case of isolated singular (i.e. non hyperbolic) point on the interior of 
one of the edges of V, call it (po,uq) = (0,0) and assume T> C {p > 0} (otherwise we apply an 
appropriate linear transformation on the conserved quantities) We investigate the propagation 
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of small nonequilibrium perturbations of the steady state of the microscopic interacting particle 
system, corresponding to the densities (po, Uq) of the conserved quantities. We prove that for a 
very rich class of systems, under proper hydrodynamic limit the propagation of these small per- 
turbations are universally driven by the system where the parameter 7 := \<& U u{po-, u o) 
(with a proper choice of space and time scale) is the only trace of the microscopic structure. 
The proof is valid for the cases with 7 > 1. 

Actually, in order to simplify some of the arguments, we impose the left-right reflection 
symmetry of the pde on the systems of interacting particles on microscopic level, see 

condition (Q in subsection 12 .21 But we note that the whole proof can be extended without this 
condition, just some arguments would be longer. 

The proof essentially relies on H-T. Yau's relative entropy method and thus, it is valid only 
in the regime of smooth solutions of the pde p.ljl . 

We should emphasize here the essential new ideas of the proof. Since we consider a low 
density limit, the distribution of particle numbers in blocks of mesoscopic size will have a 
Poissonian tail. The fluctuations of the other conserved quantity will be Gaussian, as usual. 
It follows that when controlling the fluctuations of the empirical block averages the usual large 
deviation approach would lead us to the disastrous estimate E(exp{e GAU • POI}) = 00. It 
turns out that some very special cutoff must be applied. Since the large fluctuations which 
are cut off can not be estimated by robust methods (i.e. by applying entropy inequality), only 
some cancellation due to martingales can help. This is the reason why the cutoff function must 
be chosen in a very special way, in terms of a particular Lax entropy of the Euler equation 
(|2.15|) . In this way the proof becomes a mixture (in our opinion rather interesting mixture) of 
probabilistic and pde arguments. The fine properties of the limiting pde, in particular the global 
behavior of Riemann invariants and some particular Lax entropies, play an essential role in the 
proof. The radical difference between the 7 > 1 vs. 7 < 1 cases, in particular applicability vs. 
non-applicability of the Lax-Chuey-Conley-Smoller maximum principle, manifests itself on the 
microscopic, probabilistic level. 

1.2 The structure of the paper 

In Section 12 we define the class of models to which our main theorem applies: we formulate the 
conditions to be satisfied by the interacting particle systems to be considered, we compute the 
steady state measures and the fluxes corresponding to the conserved quantities. At the end of 
this section we formulate the Eulerian hydrodynamic limit, for later reference. 

In Section 01 first we perform asymptotic analysis of the Euler equations close to the sin- 
gular point considered, then we formulate our main result, Theorem ^ and its immediate 
consequences. 

Sections El to EH are devoted to the proof of Theorem ^ 

In Section 0] we perform the necessary preliminary computations for the proof. After intro- 
ducing the minimum necessary notation we apply some standard procedures in the context of 
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relative entropy method. Empirical block averages are introduced, numerical error terms are 
separated and estimated. In this first estimates only straightforward numerical approximations 
(Taylor expansion bounds) and the most direct entropy inequality is applied. 

Section |S] is of crucial importance: here it is shown why the traditional approach of the 
relative entropy method fails to apply. Here it becomes apparent that in the fluctuation bound 
(usually referred to as large deviation estimate) instead of the tame E( exp{e GAU 2 }) we would 
run into the wild E(exp{e GAU ■ POI}) which is, of course, infinite. It is explained here what 
kind of cutoff is applied: the large fluctuations cut off can not be estimated by robust methods 
(i.e. by applying entropy inequality). Only some cancellation due to martingales can help. This 
is the reason why the cutoff function must be chosen in a very particular way, in terms of a 
particular Lax entropy of the Euler equation. The cutoff function is constructed and its key 
estimates are stated. Proofs of the lemmas formulated in this section are postponed to Section 
At the end of this section the outline of the further steps is presented. 

In SectionH3all the necessary probabilistic ingredients of the forthcoming steps are gathered. 
These are: fixed time large deviation bounds and fixed time fluctuation bounds, the time 
averaged block replacement bounds (one block estimates) and the time averaged gradient bounds 
(two block estimates). The proof of these last two rely on Varadhan's large deviation bound 
cited in that section and on some probability lemmas stated and proved in section I1UI We 
should mention here that these proofs of the one- and two block estimates, in particular the 
probability lemmas involved also contain some new (and, we hope, instructive) elements. 

Sections and El conclude the proof: the various terms arising in section [3 are estimated 
using all the tools (probabilistic and pde) developed in earlier sections. One can see that these 
estimates rely heavily on the fine properties of the Lax entropy used in the cutoff procedure. 

As we already mentioned sections OH and ^] are devoted to proofs of various lemmas stated 
in earlier parts. Section El deals with the pde estimates while Sectional is probabilistic. 

Finally in the Appendix (Section ITT|) we give some details about the pde (jl.lj) . This is 
included for sake of completeness and in order to let the reader have some more information 
about these, certainly interesting, pde-s. Strictly technically speaking this Appendix is not used 
in the proof. 

2 Microscopic models 

Our interacting particle systems to be defined in the present section model on a microscopic 
level the same deposition phenomena as the pde There will be two conserved physical 

quantities: the particle number rjj £ N and the (discrete) negative gradient of the deposition 
height Qj G Z. 

The dynamical driving mechanism is of such nature that 

(i) The deposition height growth is influenced by the local particle density. Typically: growth 
is enhanced by higher particle densities. 
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(ii) The particle motion is itself influenced by the deposition profile. Typically: particles are 
pushed in the direction of the negative gradient of the deposition height. 

The left-right reflection symmetry of the pde will be also implemented on the microscopic 
level. Actually, this is not really necessary in order to prove our main result, but without this 
assumption some of the arguments would be somewhat longer. 

2.1 State space, conserved quantities 

Throughout this paper we denote by T n the discrete tori Z/nZ, n £ N, and by T the continuous 
torus R/Z. We will denote the local spin state by we only consider the case when O is finite. 
The state space of the interacting particle system of size n is 

Configurations will be denoted 

uj_ := (ujj)j£.fn £ £l n , 

For sake of simplicity we consider discrete (integer valued) conserved quantities only. The two 
conserved quantities are 

r] : Q -> N, 

C : n -> v Z, or C : ft -> v (Z+ 1/2). (2.1) 

The trivial scaling factor vq will be conveniently chosen later (see (|2.4j) ). We also use the 
notations rjj = i](ujj), Q = ((ujj). This means that the sums ■ rjj and ^ ■ Q are conserved by 
the dynamics. We assume that the conserved quantities are different and non-trivial, i.e. the 
functions C,r] and the constant function 1 on 0, are linearly independent. 

The left-right reflection symmetry of the model is implemented by an involution 

R-.n^n, RoR = Id 

which acts on the conserved quantities as follows: 

^Rcu) = r)(u), ((Ru) = -((uj). (2.2) 

2.2 Rate functions, infinitesimal generators, stationary measures 

Consider a (fixed) probability measure tt on $7, which is invariant under the action of the 
involution R, i.e. tt{Rlo) = ir(uS). Since eventually we consider low densities of rj, in order to 
exclude trivial cases we assume that 

7r(C = 0|?7 = 0) < 1. (2.3) 
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The scaling factor vo in ()2.1|) is chosen so that 

Var(C | V = 0) = 1. (2.4) 

This choice simplifies some formulas (fixing a recurring constant to be equal to 1, see (|3.4|l ) but 
does not restrict generality. 

For later use we introduce the notation 

p* : = max{r)(w) : tt(uj) > 0}, 
u* := max{C(u;) : > 0}, 

u* := max{((w) : r]{uj) = 0, tt(uj) > 0}, 

For r, 9 G R let G(r, 0) be the moment generating function defined below: 

G{r, 9) := log ^ e Tr? ^ )+ec(w V(w). 

In thermodynamic terms G(r, 6) corresponds to the Gibbs free energy. We define the probability 
measures 

tt t j{u) := tt{u) exp(T7]{u) + 6((uj) - G{t, 9)) (2.5) 

on £1. We are going to define dynamics which conserve the quantities Y2j Vj an d Y2j Cj> posses 
no other (hidden) conserved quantities and for which the product measures 

are stationary. 

We need to separate a symmetric (reversible) part of the dynamics which will be speeded 
up sufficiently in order to enhance convergence to local equilibrium and thus help estimating 
some error term in the hydrodynamic limiting procedure. So we consider two rate functions 
r.QxQxQxQ— > R+ and s : £1 x f2 x O x — > R + , r will define the asymmetric 
component of the dynamics, while s will define the reversible component. The dynamics of the 
system consists of elementary jumps affecting nearest neighbor spins, (u}j,cjj+i) — ► {uj'-,uj'- +1 ), 
performed with rate Xr(uj, uy+i; a/, ^- +1 ) + K,s(uj, uy+i;cj'-, Wj+i)) where A, k > are speed-up 
factors, depending on the size of the system in the limiting procedure. 

We require that the rate functions r and s satisfy the following conditions. 

(A) Conservation laws: If r(ui,U2',uj' 1 , u' 2 ) > or s(ui,U2',uj' 1 ,uj2) > then 

»7(wi) + i7(w 2 ) = r?K) + ??K). 

C(^i) + CM = CK) + CK), 
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(B) Irreducibility: For every N £ [0,np*], Z £ [— nu*,nu*] the set 

[ jeT n jei n 

is an irreducible component of Q n , i.e. if lj, u/ E fi^/ ^ then there exists a series of elemen- 
tary jumps with positive rates transforming uj into u/. 

(C) Left-right symmetry: The jump rates are invariant under left-right reflection and the action 
of the involution R (jointly): 

r(Ru2, Bh)\\ Ru>' 2 , Roj'i) = r(u\, L02; uj[, u' 2 ). 
s(Ru>2, Rui'j Rw 2 , Rw'i) = s(a>i, 0J2; uj'i , u 2 ). 

(D) Stationarity of the asymmetric part: For any ui,U2,U3 G 

Q(wi, w 2 ) + Q(^2, w 3 ) + Q(w3,^i) = 0, 

where 

\ f 7T ( u; l) 7r ( u; 2) III \ I I l\\ 

(E) Reversibility of the symmetric part: For any cji , , , G ft 

ir(ui)ir(u2)s(ui,U2; u[, J 2 ) = ^(^iM^M^i; ^li ^2)- 

For a precise formulation of the infinitesimal generator on £l n we first define the map 0^ u : 
Q n -» O n for every a;', u" G fi, j G T n : 

{a;' if i = j 
c^" if i=j + l 
if + 

The infinitesimal generators defined by these rates will be denoted: 

L n m = J2 E K^,^ + i;^^")(/(0f^)-/(^)). 
K n f(ui) = E £ s(^,^+i;^^")(/(®i V ^)-/(^))- 

We denote by Xf the Markov process on the state space Vt n with infinitesimal generator G n := 
\{n)L n + n(n)K n . with speed-up factors A(n) and re(n) to be specified later 

Remarks: 
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(1) Conditions and JB| together imply that ^ • r]j and ^ • (j are indeed the only conserved 
quantities of the dynamics. 

(2) Condition (jUJ) together with Q2.2JI is implementation on a microscopic level of the left-right 
symmetry of the pde (jl.lj) . Actually, our main result, Theorem ^ is valid without this 
assumption but some of the arguments would be more technical. 

(3) Condition (|P|) implies that the product measures tt™ s are indeed stationary for the dynamics 
defined by the asymmetric rates r. This is seen by applying similar computations to those 
of CO' 12' EH1 ° r |22] ■ Mind that this is not a detailed balance condition for the rates r. 

(4) Condition (JEJ) is a straightforward detailed balance condition. It implies that the product 
measures 7r™ d are reversible for the dynamics defined by the symmetric rates s. 

We will refer to the measures 7r™ e as the canonical measures. Since YljCj an d ^jVj are 
conserved the canonical measures on £l n are not ergodic. The conditioned measures defined on 

V\z by: 

K n N , z (ui) := i$ >9 (<j\ ^ nj = N,^( j = Z,) = — 

are also stationary and due to condition {B| satisfied by the rate functions they are ergodic. 
We shall call these measures the microcanonical measures of our system. (It is easy to see that 
the measure 7r^ z does not depend on the choice of thew values of r and 6 in the previous 
definition.) 

The assumptions are by no means excessively restrictive. Here follow some concrete examples 
of interacting particle systems which belong to the class specified by conditions (|A" jl -(|E |) and 
also satisfy the further conditions (JF]), |H]), (JJ) to be formulated later. 

{—1, 0, +l}-model The model is described and analyzed in full detail in [22] and 0- The one 
spin state space is = {—1,0,-1-1} . The left-right reflection symmetry is implemented by 
R : — ► CI, Roo = -co. The dynamics consists of nearest neighbor spin exchanges and the two 
conserved quantities are r](oo) = 1 — \uo\ and ((oo) = u. The jump rates are 

r(l,-l;-l,l) =0, r(-l,l;l,-l) =2, 

r(0,-l;-l,0) =0, r(-l,0;0,-l) = 1, 

r(l, 0,0,1) =0, r(0, 1,1,0) = 1. 

and 

s(u 1 ,u 2 ;u;' 1 ,u 2 ) = 



1 if (001,002) = (oj'2,00'1) and oo\ 7^ 002 
otherwise. 

The one dimensional marginals of the stationary measures are 

1 — p ± u 

7Tp, u (0) = p, 7Tp )U (±l) = 

with the domain of variables T> = {(p, u) 6 K+ x M : p + \u\ < 1}. 
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Two-lane models The following family of examples are finite state space versions of the brick- 
layers models introduced in [21]. Let O = {0, 1, . . . , n} x {— z, — z+1, . . . , z — 1, z}, where n £ N 
and zG{|,l,|,2...}. The elements of ft will be denoted uj := Q). Naturally enough, ^ ■ rjj 
and ^ • Qj will be the conserved quantities of the dynamics. Left-right reflection symmetry is 
implemented as R : $7 — > 0, RQ) = (— c)- We allow only the following elementary changes to 
occur at neighboring sites j,j + 1: 

with appropriate rates. Beside the conditions already imposed we also assume that the one 
dimensional marginals of the steady state measures factorize as follows: 

vr(w) = 7T Q) =p(v)q(0- 

The simplest case, with n = 1 and z = 1/2, that is with 0, = {0,1} x {—1/2, +1/2}, was 
introduced and fully analyzed in [22] and |16j . For a full description (i.e. identification of 
the rates which satisfy the imposed conditions, Eulerian hydrodynamic limit, etc. see those 
papers.) It turns out that conditions (|X j) -(|E ]l impose some nontrivial combinatorial constraints 
on the rates which are satisfied by a finite parameter family of models. The number of free 
parameters increases with h and z. Since the concrete expressions of the rates are not relevant 
for our further presentation we omit the lengthy computations. 

2.3 Expectations 

Expectation, variance, covariance with respect to the measures tt^q will be denoted by E Ti #(.), 
Var T) e(.), Cov T)fl (.). 

We compute the expectations of the conserved quantities with respect to the canonical 
measures, as functions of the parameters r and 9: 

P (t,0) := E T , e (rj) = Y,v(«>>tA") = G T (r,9). 
u(t,9) := E Tj(9 (C) = C0"Kfl(w) = G 9 (r,9), 

Elementary calculations show, that the matrix-valued function 

Pt Pe \ _ ( G TT G T g 
u T u e ) \ Gg T Gee 

is equal to the covariance matrix Cov Tj g(r/, £) and therefore it is strictly positive definite. It 
follows that the function (r, 9) t— > (p(t,9),u(t,9)) is invertible. We denote the inverse function 
by (p,u) i ^ (t(p,u),9(p,u)). Denote by (p,u) i— > S(p,u) the convex conjugate (Legendre 
transform) of the strictly convex function (r, 9) *— ► G(r, 9): 

S{p,u) := sup ( P t + u9-G(t, 9)), (2.6) 

T,6 



=:G"(t,0) 
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and 

V : = {(>, it) G R+ x E : S(p, u) < 00} (2.7) 

= co{(?7,C) : ir(u) > 0}, 

where co stands for convex hull. The nondegeneracy condition (|2,3j) implies that &D f] {p = 
0} = {(0, u) : \u\ < u*}. For (/>, u) £ V we have 

r(p,u) = S p (p,u), 0(p,u) = S u (p,u). 

In probabilistic terms: u) is the rate function of joint large deviations of ■ Vji^j Cj)- in 
thermodynamic terms: <S(/9, u) corresponds to the equilibrium thermodynamic entropy. Let 

=:S"(p,u). 



l~p T'U \ I 5 pp Sp U 

^ p 9 U J V S U p S uu 



It is obvious that the matrices G"(t, 9) and S"(p, u) are strictly positive definite and are inverse 
of each other: 

G"(t, 8)S"(p, u)=I = S"(p, u)G"(t, 9), (2.8) 

where either (r, 9) = (r(p, u),9(p, u)) or (p, u) = (p(r, u(r, 9)). With slight abuse of notation 
we shall denote: 7r r{p>u)>0(piU) =: 7r PjU , =■ vr™ u , E T ( p , u ),e( P ,u) = : E p,«, etc - 

As a general convention, if f; : O m — > M is a local function then its expectation with respect 
to the canonical measure 7r™ u is denoted by 

E(p, u) := E p>u (£) = ^ . . . ,w m ,)7r PiU (a;i) • • '^(wm). 

wi,...,o; m £O m 

2.4 Fluxes 

We introduce the fluxes of the conserved quantities. The infinitesimal generators L n and K n 
act on the conserved quantities as follows: 

L n 7]i = -ip(ui,u i+1 ) +Tl>(ui- 1 ,uj i ) =: -tpi + Vi-i, 

L% = + 0(cjj_i,a;j) =: -<fo + 

K n rji= -tp s {ooi,u)i + i) + i> s {^ u uJi) =: — -01 + "01-1, 

^% = -0 s (wj,u;i + i) + ^ s (wi__i,u;j) =: -0f + 



where 



■tp{ui_,u 2 ) := K^l^M,^) (77(^2) -77(^2)) 

(£(^1,^2) := ^ r(wi,u;2;a;i,a;2)(C(w2) - C(^)) 

^(UJ!,^) ■= ^2 S ( a; l' U; 2;W 1 ,u;2)( ? ?(^2) - 7 ?(^2)) 

<j) s {u)l,U2) := ^2 sO^Ij^MMHCO^) - C(^2)) 
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(2.c 



(2.10) 



Note that due to the left-right symmetry and conservations, i.e. (|2.2|) and conditions and 
microscopic fluxes have the following symmetries: 



^(o>i,o> 2 ) = (j)(Ru2,Rui), 

■0(o>i,O> 2 ) = -lf)(Ru}2,RuJ\). 

In order to simplify some of our further arguments (in particular, see (|7.7|) in subsection 
17. 3j) we impose one more microscopic condition 

(F) Gradient condition on symmetric fluxes: The microscopic fluxes of the symmetric part, 
defined in (|2.1U|) satisfy the following gradient conditions 

■^ s (oi,o 2 ) = k(u}i) - k(uj 2 ) =■ - k 2 ( 2 -ll) 

S (oi,O 2 ) = X(^l) - X{U2) =■ Xl - X2- 

Remark: (1) This is a technical assumption (referring actually to the measure tt) which sim- 
plifies considerably the arguments of subsection 19.21 The symmetric part K n has the role of 
enhancing convergence to local equilibrium. Its effect is not seen in the limit, so in principle 
we can choose it conveniently. Without this assumption we would be forced to use all the non- 
gradient technology developed in [2^] (see also JOJ), which would make the paper even longer. 
(2) It is easy to see that 7/(o>i) = 7/(0)2) = implies t/> s (o>i, 02) = and thus (by choosing a 
suitable additive constant) o> 1— > k(uj) can be chosen so that 

ri(u) = => k(u) = 0. (2.12) 

The macroscopic fluxes are: 
tt(p,ti) := E p , n (V0 



r(oi,o; 2 ;o' 1 ,o; 2 ) (r/(a; 2 ) - 77(0-2 ))ir PjU (ui) ^(02), 



^1,^2, 



$(p,u) := E AU (0) 



(2.13) 



^2 r(oi,02;o' 1 ,o 2 )(C(^2) ~ C (^2 ))?!>,« (wi)tt p ,u (^2) • 



These are smooth regular functions of the variables (p, u) £ T>. Note that due to reversibility 
of K n , for any value of p and u 

^ P , u m = = E Piti ((/> s ). 

(These identities hold true without assuming condition 
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For later use we mention here that according to |22| . the macroscopic fluxes ^(p, u) and 
$(p, u) satisfy the following Onsager reciprocity relation 

tf u (p, u)Var PjU (C) - $ u (p, u)Cov p>u (?7, C) = (2.14) 

*p(p, «)Var P)U (r?) - * p (p, u)Cov p>u (r/, £). 

For the concrete examples presented at the end of subsection 12,21 the following domains T> 
and macroscopic rates are gotten: 

{-1,0,+!} -model : 

V = {(p, u) £ M+ x M : p + |u| < 1} 
^(p, = pu 
$(p,u) = p + u 2 . 

Two /ane models with n = 1: 

V = {(p,u) e R + X R : p < 1, |u| < z} 

V(p,u)=p(l-p)i/)(u) 

$(p,u) = ip (u) +ptp 1 (u), 

where is °dd, while ipo(u) and <pi(it) are even functions of u, determined by the jump 
rates of the model. In the simplest particular case with z = 1/2 

*(p,ii) = p(l - p)u 
$(p,«) = (p- 7 )(l-u 2 ), 

where 7 6 R is the only model dependent parameter which appears in the macroscopic fluxes. 
For details see [22"] . 

2.5 The hdl under Eulerian scaling 

Given a system of interacting particles as defined in the previous subsections, by applying Yau's 
relative entropy method (see [22] or the monograph JOJ ) , one shows that under Eulerian scaling 
the local densities of the conserved quantities p(t,x), u(t,x) evolve according to the system of 
partial differential equations: 

dtp + d x *(p,u) = 

2.15 

d t u + d x <f>(p,u) = 

where ^(p,u) and $(p, u) are the macroscopic fluxes defined in (|2.13|) . 

The precise statement of the hydrodynamical limit is as follows: Consider a microscopic 
system which satisfies conditions of subsection l2.2l Note that condition (JFJ of subsection 
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12.41 is not assumed. Let u) and u) be the macroscopic fluxes computed for this system 
and p(t,x),u(t,x) x G T, t £ [0, T] be smooth solution of the pde (|2.15j) . Let the microscopic 
system of size n be driven by the infinitesimal generator 

G n = nL n + n 1+s K n , 

where 5 € [0, 1) is fixed. This means that the main, asymmetric part of the generator is speeded 
up by n and the additional symmetric part by n 1+<5 . Let pfi be a probability distribution on 
£l n which is the initial distribution of the microscopic system of size n, and 

the distribution of the system at (macroscopic) time t. The local equilibrium measure (itself 
a probability measure on f2 n ) is defined by 

:= II %(fi),«i(fi)- 

This measure mimics on a microscopic scale the macroscopic evolution driven by the pde 1)2.15(1 . 

We denote by H(p™\7r n ), respectively, by H{p^\v^) the relative entropy of the measure 
pi with respect to the absolute reference measure 7r n , respectively, with respect to the local 
equilibrium measure v%. 

The precise statement of the Eulerian hydrodynamic limit is the following 

Theorem. Assume conditions Q) )-fS )) and let 5 £ [0, 1) be fixed. If 

H(p%\vZ)=o(n) 

then 

H(p?\v?) = o(n) 

uniformly for t £ [0, T] . 

Remark: Note that due to finiteness of the state space O the condition 
holds automatically. 

The Theorem follows from direct application of Yau's relative entropy method. For the proof 
and its direct consequences see |2fi| . |l()j or [221- For the main consequences of this Theorem 
see e.g. Corollary 1 of |2"2"] . 
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3 Low density asymptotics and the main result: hydrodynamic 
limit under intermediate scaling 

3.1 General properties and low density asymptotics of the macroscopic fluxes 

The fluxes in the Euler equation l|2.15|) are regular smooth functions of in (p, u) G T>. 
From the left-right symmetry of the microscopic models it follows that 

$(p, -u) = $(p, u), V(p,-u) = -*(p,u). (3.1) 

It is also obvious that for u E [— u*, it*] 

tf(0,it)=0. (3.2) 

We make two assumptions about the low density asymptotics of the macroscopic fluxes. 
Here is the first one: 

(G) We assume that fy pu (0, 0) ^ 0. Actually, by possibly redefining the time scale and orien- 
tation of space, without loss of generality we assume 

¥^(0,0) = 1. (3.3) 



From the Onsager relation 1)2.14(1 and obvious parity considerations it also follows that 

$ p (0, 0) = ¥ pu (0, 0)Var 0i0 (C) = 1. (3.4) 



Note, that here we rely on the choice 1)2.4(1 of the scaling factor vq in (|2.1|) . 
We denote 

7:=^*««(0,0). (3.5) 

Our results will hold for 7 > 1 only. 
From (|3~Tj) and it follows that 

$ u (0, u) - tt p (0, u) = (2 7 - l)u + 0(\u\ 3 ). (3.6) 

The second condition imposed on the low density asymptotics of the macroscopic fluxes is: 

(H) For u E u^O 

*„(0,tt) -%(0,u) #0, (3.7) 

$ p (0,u)^0, * pu (0,u)^0 (3.8) 
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Remarks: (1) (JHJ is a very natural nondegeneracy condition: if ^ ptt (0, 0) vanished then in the 
perturbation calculus to be performed, higher order terms would be dominant and a different 
scaling limit should be taken. 

(2) Due to (|3.1j) . ()3.3j) and (|3.6j) conditions ()3.7|) . (|3.8j) hold anyway in a neighborhood of u = 0, 
and this would suffice, but the forthcoming arguments, in particular the proof of Lemmas [2] 
and |31 would be less transparent. We assume condition JHJ for technical convenience only. 
Condition 1)3. 7 j) amounts to forbidding other non-hyperbolic pointson dT> n {p = 0}, beside the 
point (p,u) = (0,0). Condition ()3.8j) reflects the natural monotonicity requirements (i) and (ii) 
formulated about the microscopic models at the beginning of Section [2j 

We are interested in the behavior of the pde near the isolated non-hyperbolic point (p, u) = 
(0,0). The asymptotic expansion for p + u 2 <C 1 of the macroscopic fluxes and their first partial 
derivatives is 

u) =pu(l + 0{p + u 2 )) , $(p, u) = (p + 7 u 2 ) (1 + 0(p + u 2 )) , 

* p (p,u) = u(l + 0(p + u 2 )), <5> p (p,u) = l + 0(p + u 2 ), (3.9) 

* u (p,u) = p{l + 0{p + u 2 )), $ u (p,u) = 2 1U (l + 0(p + u 2 )). 

We are looking for "small solutions" of the pde (|2.15|) : Let po(x) and uq(x) be given profiles 
and assume that p £ (t,x), u £ (t,x) is solution of the pde (|2. 15|) with initial condition 

p e (0, x) = e 2 po(x), u £ (0,x) = euq(x). 

Then, at least formally, 

e~ 2 p e (e~ 1 t,x) — > p(t,x), e~ 1 u e (e~ 1 t,x) — > u(t,x), 
where p(t,x), u(t,x) is solution of the pde ([l.lj) with initial condition 

p(0, x) = p (x), u(0, x) = uq{x). 

3.2 The main result 

The asymptotic computations of subsection l3.1l suggest the scaling under which we should derive 
the pde as hydrodynamic limit: fix a (small) positive (3 and choose the scaling 

MICRO MACRO 

space nx x 

time n 1+/3 t t 

particle density n~ 2 ^p p 

'slope of the wall' n~^u u 

Ideally the result should be valid for < (3 < 1/2 but we are able to prove much less than that. 

Choose a model satisfying the conditions (0)-(0 of section El and conditions (jOllrlj) of 
subsection 13.11 and let 7 be given by (|3.5|) . corresponding to the microscopic system chosen. 
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Let the microscopic system of size n (defined on the discrete torus T n ) evolve (on macroscopic 
time scale) according to the infinitesimal generator 

G n = n 1+ ^L n + n 1+ P +5 K n . 

with /3 > and some further conditions to be imposed on j3 and 5 (see Theorem Denote by 
the true distribution of the microscopic system at macroscopic time t: 

where p,Q is the initial distribution. 

We use the translation invariant product measure 

71 - =7r n-^,0 

as absolute reference measure. Global entropy will be considered relative to this measure, 
Radon-Nikodym derivatives of and the local equilibrium measure to be defined below, 
with respect to 7r n will be used. 

Given a smooth solution [p(t, x), u(t, x)) , (t, x) £ [0, T] x T, of the pde define the local 
equilibrium measure on Q n as follows 

v t '■= II K n -2f) p( t £) in -/3 u(t !)■ (3.10) 

This time-dependent measure mimics on a microscopic level the macroscopic evolution governed 
by the pde (fT~l"i . 

Our main result is the following 

Theorem 1. Assume that the microscopic system of interacting particles satisfies conditions 
QP-~EP of subsubsections \2. 6 A ~~~~| and the uniform log-Sobolev condition |7|) of subsection \6.& 
Additionally, assume that the macroscopic fluxes satisfy conditions f?7j) . flT)) of subsection \H.f\ 
and 7 > 1. Choose (3 £ (0, 1/2) and 5 £ (1/2, 1) so that 

25- 8/3 >1 and 5 + 3p<l. (3.11) 

Let (p(t, x), u(t, x)) , (t,x) £ [0, T] x T, 6e smooth solution of the pde M.l\) . such that 

inf^gf p(0, x) > and let , t £ [0, T] be the corresponding local equilibrium measure defined 
in 137771) . 

Under these conditions, if 

H{tf | iff) = o(n 1 ^) (3.12) 

then 

H(tf\v?)=o(n 1 - 2 P) (3.13) 

uniformly for t £ [0, T] . 
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Remarks: 

(i) From (|3. 12|) via the identity (|4,5|) and the entropy inequality it also follows that 

H{^\ii n ) =0{n 1 - 2 < 3 ). (3.14) 



See the beginning of subsection 14.21 

(i) If 7 > 3/4, in smooth solutions vacuum does not appear. That is inf x . e T p(0, x) > implies 
itlf ((,i)e[0,T]xTP(^) > 0. 

(ii) Although for the {—1,0, +l}-model we have 7 = 1, our proof can also be extended to cover 
this model. Actually, in that case the proof is much simpler, since the Eulerian pde is equal to 
the limit pde (|l.lj) and thus the cutoff function (see Section EJ) can be determined explicitly. 

Corollary 1. Assume the conditions of Theorem^ Let g : T — > E be a test function. Then 
for any t € [0, T] 

(i) 



1 Y] g(-)Vj(t) ^ g(x)p(t,x)dx, 
n It 

TA'W * / g(x)u(t,x)dx. 



/? 

jeT 
(ii) 

H{^\n n )-H{^\n n )=o(n 1 ^). 
See the proof of Corollary 1 in [221 • 

4 Notations and general preparatory computations 

This section completely standard in the context of the relative entropy method. So we shall be 
sketchy. 

4.1 Notation 

We denote 

: = n-^(H(^\-K n )-H(^\n n )). 
We know a priori that 1 1— > s n (t) is monotone increasing and due to ()3.14|) 

s n (t) = 0(1), uniformly for t G [0, oo). (4.1) 
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In fact, from Theorem^it follows (see Corollary^) that as long as the solution p(t,x),u(t,x) 
of the pde is smooth 

s n (t) = o(l), uniformly for t G [0, T). 

For (p,u) G (0,oo) x (—00,00) denote 

T n (p, u) := r(n" 2 V n~ p u) - r(n" 2/3 , 0) 

n (p,u) :=n^(n- 2/3 />,n~^)- 

Note that, for symmetry reasons 9(n~ 2 ^,0) = 0. Mind that r is chemical potential rather than 
fugacity and for small densities the fugacity A := e T scales like p, i.e. r(n _2 ^,0) ~ —2(3 log n. 
If p > and u G R are fixed then T n (/3, it) and # n (/0, u) stay of order 1, as n — > 00. 

Given the smooth solution p(t,x),u(t,x), with p(t,x) > we shall use the notation 

T n (i,x) :=r n (/3(i,x),n(t,2;)), 

r(t,x) :=r(p(t,x),n(t,x)), 

v(t, a;) := log /?(£, x). 

The following asymptotics hold uniformly in (t, x) G [0, T] x T: 

r n (t,x)= v(t,x) + 0(n~ 2 P), 6 n (t,x)= u(t,x) + 0{n- 2 < 3 ) 
d x r n (t, x) = d x v(t, x) + 0(n~ 2 P), d x 9 n (t, x) = d x u(t, x) + 0(n~ 2fi ) (4.2) 
8 t r n (t, x) = dt v(t, x) + 0(n~ 2 P), d t 6 n (t, x) = dt u(t, x) + 0(n- 2 P) 

The logarithm of the Radom-Nikodym derivative of the time dependent reference measure 
with respect to the absolute referencee measure 7r n is denoted by 

= E{ rn (*^H+^ r (*^)o ( 4 - 3 ) 

-G(r n (t, 3 -) + r(n- 2 ^, 0), n~H n (t, J -)) + G(r(n" 2/3 , 0), 0)| 
n n > 

4.2 Preparatory computations 

In order to obtain the main estimate (|3.13|) our aim is to get a Gronwall type inequality: we 
will prove that for every t G [0, T] 

h n (t) - h n (0) = f d t h n {s)ds <C f h n (s)ds + o(l), (4.4) 
Jo Jo 

where the error term is uniform in t G [0, T]. Because it is assumed that h n (0) = o(l), the 
Theorem follows. 
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We start with the identity 

H{£\u?) - H(n?\ir n ) = - [ tfdtf. (4.5) 

From this identity, the explicit form of the Radon-Nikodym derivative (|4.3j) . the asymptotics 
(|4.2j) . via the entropy inequality and (|3.12|l the a priori entropy bound (|3. 14f) follows indeed, as 
remarked after the formulation of Theorem Q 
Next we differentiate (j4.5|) to obtain 

d t h n (t) = - [ (n zp L n ff + n 3l3+s K n f? + rT 1+2P d t ff) d$ - d t s n {t). 

(4.6) 

Usually, an adjoint version of Q4.6|) is being used in form of an inequality. In our case this 
form is needed. We emphasize that the term —dts n {t) on the right hand side will be of crucial 
importance. 

We compute the three terms under the integral. 

i 3 ^Ty = - E w*, 3 -> z ^3 + - E d * u ^ L ) nW ^ ( 4 -7) 

n n n ^-^ n 

+A\ (t, uf) + A%{t,uj_) + A% (i, u) +A2(t,u), 
where Af(t,uj_), i = 1, . . . ,4 are error terms which will be easy to estimate: 

n ^-^ n n 

A n 2 {t^) := -V(8/(^)-H-))A. 

n ' n n 

jr'GT™ 

A$(t,u) := iV(VT(^)-9/(^))^, 
n n n 

Here and in the sequel V™ denotes the discrete gradient: 

V n f(x) :=n(f(x + l/n)-f(x)). 

See subsection 14.41 for the estimate of the error terms A™(t,uj), j = 1, . . . , 12. 

Next, 

n 3f3+s K n f?(u/) = n - 1+3l3+s - V ( (V n ) 2 r n {t, ?-) Kj + (V n ) 2 d n (t, i-\ x .) 

= : A^u) (4.8) 
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is itself a numerical error term. Finally 

n-i+Wdtfn^ = I V (n 2 %--p(t,^)) (4.9) 



n ^— ' n v n 

+Ag(t,u;) + ,4?(i,u/), 



where 



^(^) == -E (^ n (*.-)-^,-))(n 2 %--p(t,^-)), 
n n n n 

A^) := \ E (d t e n (t, J -) - d t u(t, 3 ~)) (n% - u(t, 3 -)). 

ft IV IV IV 

are again easy-to-estimate error terms. 
4.3 Blocks 

We fix once and for all a weight function a : K — > R. It is assumed that: 

(1) a{x) > for x 6 (—1,1) and a(x) = otherwise, 

(2) it has total weight J a(x) dx = 1, 

(3) it is even: a(— x) = a(x), and 

(4) it is twice continuously differ entiable. 

We choose a mesoscopic block size I = l(n) such that 

1 < n (l+<5+5/3)/3 ^ /( n ) ^ n 5-/3 ^ n _ ( 410 ) 

This can be done due to condition (|3.11|) imposed on j3 and 5. 
Given a local variable (depending on m consecutive spins) 

& = £i(<^) = £( w i> • • • > ^i+m-l), 

its Woc/c average at macroscopic space coordinate x is defined as 



i ^ v * 

Since / = Z(n), we do not denote explicitly dependence of the block average on the mesoscopic 
block size I. 

Note that x i— ► £ n (x) is smooth 



j 
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and it is straightforward that 



sup sup 



d x C(u, x) < C [ max £(u>i, . . . , u Vv/ 



n 



(4.12) 



For a more sophisticated bound on 
We shall use the handy (but slig 



d x l n {u, x) 



see (jfrm 



rtly abused) notation 

?{t,x) :=f>(*M- 

This is the empirical block average process of the local observable £j. 

For the scaled block average of the two conserved quantities we shall also use the notation 

fT{t,x) ■.= n 2P rf l {t,x), u n (t,x) :=n^C{t,x). 

Introducing block averages the main terms on the right hand side of (|4,7|) and Q4.9|) become: 



- V d x v(t, l W% + - E M*> L > w ^ 

n ^— ' n n n 



v d x v(t, 3 -)n^r{t) + 1 v axt, -)™ 2/3 ^ n (-) 

^ — * 77, r) r) ^ — * 

j£T n 



n 



n n 



n 



jeT n 

+A%(t,u) + AUt,UL), 



n 



respectively 

- d t v(t, ^)(n 2 %- - p(t, + dtu(t, J -)(n% ~ u(t, ?-)) 



n 



n 



n n 



n 



1 



n ^-^ n n n ' 



+AUt,uL) + A^(t,uD, 



where the error terms are 
A^ui) 



1 

n 



i - ^ 



71' 



J? 



E(M^)-jE 



(4.13) 



i E M<--)W-)-pft-)) ( 4 - 14 ) 
n n n n 
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These error terms will be estimated in subsection 14.41 

Since [0, T] x T 3 (t,x) i— > (p(t,x),u(t,x)), is a smooth solution of the pde we have 

3tv = —ud x v — d x u, dtu = —pd x v — jud x u. 

Inserting these expressions into the main terms of (|4,14|) eventually we obtain for the integrand 
in l|4~H|) 

W L n f t l {ui) + n 3/3+<5 iTY"(w) + n- 1+2/3 dtf?(ui) = (4.15) 



n 



I d x v(t, l){n^r( 3 -) - p(t, 3 -)u{t, I) 
n n V n n n 

- u(t, 3 -){n 2 ^( L ) - P(t, ^)) - p(t, ^)(n^ n (-) - <t, 

n n n n n n ' > 

+ - V d x u(t, t)f n ^(l) - ( p (t, I) + 7 u(i, 3 -f) 
n n I n n n 

- (" 2/ V (-) ~ Pit, -)) - 2 7 u(t, J -){n^C( J -) - u(t, 3 - 
v n n ' n n n 

12 



k=l 



where 



A n® : = ~ E {{d x v)pu+{d x u){p + 1 u 2 )){t, 3 -) 
4.4 The error terms A%, k = 1, ... ,12 

Lemma 1. There exists a finite constant C, such that for any n £ N, t 6 [0, T] and /or any 
sequence of real numbers bj, j = 1, . . . , n i/ie following bounds hold: 

Mi E m,) * c "" e E <>> + ,/^E"^ 



n * — ' / \ n 



«W (; E Wa) £ ^ f^E»i + E «?) 

*v? (I E w>) ^ ( ~ E h + E q J ( 4 - 19 ) 
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Proof. The proof relies on the entropy inequality 



(4.20) 



-^if (/4V n ) + ^ log ( exp {7 E 6 i fe " (6)) }) , 

II A rlTn 



where £j stands for either of 77^ , (j, i/ij or We note that all these variables are bounded and 
I E^n (r]j) I < Cn" 2/3 , Var^n < Cn~ 2 ^, 

I e tt« (Cj) I = 0, Var,n(C J )<C, 
I E T n (^) I = 0, Vaivn < Cra" 2 ^, 

I E^n (^) I < C, Vaivn < C. 

From these bounds and the entropy inequality (|4.2U[) the statement of the lemma follows directly. 

□ 

Now we turn to the estimates on the error terms. We use the bounds ()4.16j) . (|4.17|) . (|4.18f) 
and (|4.19|) of Lemma ^ the asymptotics (j4.2|) and uniform approximation of d x of smooth 
functions by their discrete derivative V n . Straightforward computations yield 



E 



ft 



E 



ft 



E 



ft 



E 



ft 



E 



ft 



E 



ft 



• < C(n~ 1+/3 + n~P) = o(l), 
<C(n~ 1+2 ' 3 + n-' 3 ) = o(l), 

.(Am) <Cn- 1+ P = o(l), 
\Al{t)) < C(n- 1+2/3 + n- 1+/3 ) =o(l), 
■ (A2(t)) < C(n~ 1+ P +S + n~ 1+2 ' 3 + s ) = o(l), 
\Alit)) < Cn-W = o(l), 

(A?(i)) < C7i~ 2/3 = O(l), 

• (^s (*)) < C(n" 1+/3 + r/r 1 + n" 1+/3 Z) = o(l), 
. (^(t)) < C(n" 1+2 ^ + n^r 1 + n" 1+/3 /) = o(l), 

(^io(*)) < C ( n " X + rl + n ~^) = ol 1 )' 

(A^it)) <c(r l + n~H) =o(l). 

Finally, A™ 2 (f) ls a simple numerical error term (no probability involved): 

A? 2 (t) < CrT 1 = o(l). 



E 



ft 



E, 



E 



ft 



E 



ft 



E 



ft 
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4.5 Sumup 

Thus, integrating (|4.fi|) . using (|4.15j) and the bounds of subsection 14.41 we obtain 

h n {t)= [ A n (s)ds+ I B n (s)ds- s n (i) + o(l), (4.21) 
Jo Jo 

where 

A n (t) := E(i Yl {{d x v){n^r - pu - u(fT - p) - p(u n - u)}}(t, J -)) (4.22) 

jr'GT" 

and 



B n (t) := E(i Yl { (MK^ - (p + ju 2 ) - (p» - p) - 2 7 n(u" - «)}}(*, ^)) (4.23) 

The main difficulty is caused by .A n (i). The term B n (t) is estimated exactly as it is done in 
|21j for the one-component systems: since Q(p,u) = p + 7 u 2 is linear in p and quadratic in u 
no problem is caused by the low particle density. By repeating the arguments of we obtain 

f B n (s)ds<C [ h n (s)ds + o(l). (4.24) 
Jo Jo 

In the rest of the proof we concentrate on the essentially difficult term A n (t) . 

5 Cutoff 

We define the rescaled macroscopic fluxes 

* n {p, u) := n 3f3 ^{n- 2f3 p, rC P u), <f> n (p, u) := n 2/3 $(n~ 2/ V, n^u). (5.1) 
defined on the scaled domain 

V n := {{p, u) : (n- 2f3 p, rT P u) G V}. (5.2) 

The first partial derivatives of the scaled fluxes are 

^(p, u) = n 13 ^ p (n~ 213 p, n'^u), &p(p, u) = ^ p (n~ 2l3 p, n'^u), 

*2(/o, u) = n 2 P$> u (n- 2 Pp, n~P U ), $£(p, u) = ji^^i" 2 ^, rT fi u). 
For any (p, u) G E + x E 

lim ^ n (p,u) = pu, lim & n (p, u) = p + ju 2 , 

n— >oo n— »oo 

lim ^™(p,u) = u, lim $"(p,u) = 1, (5.4) 

n— >oo ' n— >oo r x ' 

lim = lim K(P, U ) = 2 JU. 

n— >oo n— >oo 

The convergence is uniform in compact subsets of M4. x R 
Note that 

s), S»(t, x)) = n 3 ^(^(t, x), C" (t, x)), 
^(i, x), u"(t, x)) = n 2 ^(fT(t, x), C(i, x)). 



(5.3) 
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5.1 The direct approach — why it fails? 

The most natural thing is to write the summand in A n (t) as 

n 3 ^ n -pu- u(pr - p) - p(u n -u)= (5.5) 

By applying Varadhan's "one block estimate" and controlling the error terms in the Taylor 
expansion of the first two terms on the right hand side can be dealt with. However, the 
last term causes serious problems: with proper normalization, it is distributed with respect to 
the local equilibrium measure f", like a product of independent Poisson and Gaussian random 
variables, and thus it does not have a finite exponential moment. Since the robust estimates 
heavily rely on the entropy inequality where the finite exponential moment is needed, we have 
to choose another approach for estimating A n (t). 

Instead of writing plainly ()5.5j) . we introduce a cutoff. We let 

M > sup{p(t,x) V \u(t,x)\ : (t,x) E [0, T] x T}. 

The value of M will be specified by the large deviation bounds given in Proposition |2] (via 
Lemma EJ). 

Let I n ,J n : R + x M —> R be bounded functions so that 

jn + j" = l5 

I n (p,u) = l for p\/\u\<M, 
I n (p,u) = for 'large' {p,u). 

The last property will be specified later. 

We split the right hand side of (|5.5|) in a most natural way, according to this cutoff: 

n 3 ^ n -pu- u(pr -p)- p(u n -u)= (5.6) 
n 3/3 ^ n J n (p n , u n ) - {p^u + pu n - pu)J n {p^,u n ) 

+(^ n (p n ,n") - p^u n )I n {pr,u n ) + (p n - p){u n - u)I n {p^,u n ) 

The second term on the right hand side is linear in the block averages, so it does not cause any 
problem. The third term is estimated by use of Varadhan's one block estimate. The fourth term 
is Taylor approximation. Finally, the last term can be handled with the entropy inequality if 
the cutoff I n (p,u) is strong enough to tame the tail of the Gaussian x Poisson random variable. 

The main difficulty is caused by the first term on the right hand side. This term certainly 
can not be estimated with the robust method, i.e. with entropy inequality: we would run into 
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(5.8) 



the same problem we wanted to overcome by introducing the cutoff. The only way this term 
may be small is by some cancellation. It turns out that the desired cancellations indeed occur 
(in form of a martingale appearing in the space-time average) if and only if 

J n (p,u)=S;(p,u), (5.7) 

where S n (p,u) is a particular Lax entropy of the scaled Euler equation 

r d tP + d x ^ n (p,u) = o 
\ d tU + d x ^ n ( P ,u) = o, 

with ^ n (p, u) and $ n (p, u) defined in ()5.1[) . That is S n is solution of the pde 

Ks n pp + {K - n)s? u - $> n P s n uu = 0- (5.9) 

5.2 The cutoff function 

In the present subsection we construct the cutoff function (|5.7j) and we state some estimates 
related to it. These bounds will be of paramount importance in our further proof. They are 
stated in the technical Lemmas 12 El and 0] The proof of these lemmas is pure classical pde 
theory and it is postponed to section 

First, in subsubsection IB~ 2 . 1 1 we formulate our construction and estimates in terms of Lax 
entropies of the unsealed Euler equation (|2.15j) . Then in subsubsection 15.2.21 we rescale these 
estimates in order to get the necessary bounds on S n and its derivatives. 

5.2.1 

A Lax entropy/flux pair S(p,u), F(p,u) of the system ()2.15|) is solution of the system of pdes 

F p = * P S P + $ P S U , F u = ^ U S P + § U S U , (5.10) 

defined on T>. In particular the Lax entropy S(p, u) solves the pde: 

V U S PP + {$u-* P )S fm -$ p S uu = 0, (5.11) 

The linear pde (|5.11j) is hyperbolic in T>. One family of its characteristic curves are solutions 
of the following ODE, meant in the domain T>: 



dp _ V - *") 2 + 4$ptt u-($u-*p) 

, (5.1ZJ 



du 2®p 



The other family is obtained by reflecting u to —u. 

First we conclude that the line segment T>f]{u = 0} is not characteristic for the hyperbolic 
pde (j5.11|) . That is: is intersects transversally the characteristic lines defined by the differential 
equation ()5.12() . Indeed, from the Onsager relation (|2,14[) and obvious parity considerations it 
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^ 1 /2 

follows, that the right hand side of (j5.12j) restricted to {u = 0} becomes (Var r) o(77)/Var ri o(C)) 
and this expression is obviously finite for r G [0, p*). It follows that the Cauchy problem (|5.11|) . 
with the following initial condition: 

S(r, 0) = s(r), S u (r, 0) = 0, r G [0, p*) (5.13) 

is well posed. 

In our concrete problem the function s(r) will be chosen as follows: we fix < r < r, and 
define 

'0 if r G [0,r), 

r log (r/r) — (r — r) 



s(r) 



log (r/r) 
r - r 

it r G [r, oo 



ifrG[r,r), (514) 



log (r/r) 

Note that s(r) and s'(r) are continuous. 

We first analyze the global structure of the characteristic curves. Due to the assumption 
(|H|) imposed on-, and regularity of the flux functions $ and there exists some po > such 
that the ODE (f5~T2^ is regular in {(p, u) G V : p < p and (p,u) / (0,0)}. We shall not be 
concerned about what happens outside this strip. Denote by o~(u; r) the solution of the ODE 
(|5,12[) with initial condition cr(0; r) = r. 

Lemma 2. There exist constants < C 2 < C\ < 00 and ro > such that for any r G [0, ro] 

T + C\y/ru < a(u; r) < r + C 2 \fru if u < 0, 

/ 47-3 1 \ (5.15) 

r < a(n; r) < r + C± I \fru A r 4 t~ 2 u 2 ^- 1 J if u > 0. 

The inequalities are valid as long as (a(u;r),u) G "D. T/ie map u 1— > <r(u;r) is regular and 
monotone increasing. 

See subsection 19.11 for the proof of this lemma. 

For r < ro we partition the domain T> in three parts as follows 

Vx(r) := {{p,u)€V:p<a{-\u\;r)} 

V 2 (r) := {{p,u) £P : p> a{\u\;r)} 

V 3 (r) := {{p,u)eV:a(-\u\;r)<p<a(\u\;r)} 

= V\(V 1 (r)UV 2 (r)). 

See Figure 1 for a sketch of the domains 2?i(r), T> 2 (r), T>^{r). 
From Lemma [21 it follows that 

{{p,u) : < p < r - C X y/r\u\} C X>i(r) C {(p,u) : < p < r - C 2 Vr\u\}, 

(5 16) 

{(p,u) : r + C ly fr\u\ < p} C £> 2 (r) C {(p,u) : r < p}, 
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and for < r < r' < ro 



V x (r) C V 1 (r'), V 2 {r') C V 2 (r) 



(5.17) 




U 

-u* U* 

Figure 1:V 1 ,V 2 ,V 3 

From now ro is fixed for ever and we denote 

V:=Vx(ro). 

This domain is a rectangle in characteristic coordinates with diagonal Df){u = 0}, as opposed 
to T> which may not be a full characteristic rectangle. (Actually, choosing the characteristic 
coordinates in a natural symmetric way, z(p,u) = w(p, —u), the domain T> is a square in char- 
acteristic coordinates.) Note that V^^q) n {u > 0} and X^ro) D {u < 0} are also characteristic 
rectangles. 

Next we turn to the construction of a particular family of Lax entropies which will serve for 
obtaining the cutoff functions needed. We fix < r < r < tq. and define S : T> — > R as follows: 



(i) In T>: S(p,u) is solution of the Cauchy problem (|5.11j) + ()5.1Hj) with s(r) given in (|5.14j) . 
Note that 

' if (p,u) e Vi{r) C V, 



S(p,u) 



P 



(ii) InD 2 (f): 



S(p,u) := p 



log (r/r) 



r — r 



if ( P ,u) ev 2 (r)nv. 



if ( P ,u)ev 2 ( 



log(r/r) 

Note that there is no contradiction: in D D T> 2 (r), (JIJ) yields the same expression. 



(5.19) 
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(iii) In T> 3 (r) \ D: S(p,u) is denned as solution of the Goursat problem (|5.11|) with boundary 
conditions on the characteristic lines dVr\V 3 (r), respectively, dT>2(r) \ V provided by JlJ), 
respectively, ffTT|) . 



Note that S(p, u) is solution of the pde ()5.11j) . globally in T>. 
We denote 

V 3 (r,r) := V \ (l»i(r) U V 2 (r)) 

= (2?i(r) n V 2 {r)) U (2?i(r) n V 3 (r)) 

u(V 3 (r)nV 2 {r)) U (V 3 (r) n V 3 (r)) . 



P 




u 



Figure 2: ^(r, ^) 

The following lemma provides the necessary bounds on the partial derivatives of S(p, u) (up 
to second order) in the domain T> 3 (r, r). 

Lemma 3. There exists a constant C < oo such that for any < r < r < ro and S(p, u) 
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defined as above the following global bounds hold: 



\S p (p,u) - l V2{7) (p,u) 
\S u (p,u) 

\S PP (p,u) 

\S pu (p,u) 

\S UU (P,U) 



< CTL V3(z ^(p,u 
C(Vf-,/r 



< 



< 



< 



< 



log (r/r) 
C 1 



log(r/r) r + p 
C 1 



log(r/r) + a/p + 1^ 

C _ , , 
log (r/r J v y 



lr> 3 (r,r)(P, «), 

lx> 3 (z:,r)(P> M ): 



(5.20) 
(5.21) 

(5.22) 

(5.23) 

(5.24) 



This lemma is proved in subsection 19.21 

Beside the bounds on the partial derivatives of S(p,u) we shall also need a bound on the 
function F(p,u) — ^f(p,u)S p (p,u). From ()5.1Uj) and (|5.18|) it follows that 



F(p,u) 



f if (p,u) £Di(r), 
I if (p,u) £V 2 (f). 

Lemma 4. With the assumptions and notations of Lemma\^ 

\F(p,u) - ^{p,u)Sp(p,u)\ < - — -——(r + u 2 )l v irT) (p, 

log(r/r) v - ' 

See subsection 19.31 for the proof of this lemma. 



u 



(5.25) 



5.2.2 

The scaled functions S n (p,u), F n (p,u) are defined on the scaled domain T> n given in 1)5 .2 J) , as 
follows: fix0<r<r<oo and define the unsealed Lax entropy/flux pair as in the previous 
section but with downscaled initial conditions 

S{r,0) =n- 2/3 s(n 2l3 r), S u (r,0)=0. r€[0,p*), (5.26) 

with the function r i— > s(r) given in (|5.14jl . Now, define the pair of scaled functions S n ,F n : 
V n -> R as 

S n (/?, u) := n 2/3 S(n- 2/ V, F n (p, u) := n 3/3 F(n~ 2/3 p, n^u). (5.27) 
It is straightforward to check that S n , F n form a Lax entropy/flux pair of the pde (|5.8|) : 

F; = *^ + $^, F- = ^p+«, (5-28) 
in particular S n solves the pde ()5.9|) . 
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We partition the scaled domain 

V n = 2?J(r) U V%(r) U V n A (r,r) 

with the partition elements 

E>«(r) := {(p,it) G X>„ : («" 2 Vn^«) G Pi(n-^r)} 

V%(r) := {(p,u) G P n : {rC^p^u) G P 2 (n" 2/3 r)} 

P"(r,r) := {(p,u) G £>„ : (n~ 2f3 p, rT p u) G V 3 (n~ 2 ^r, n~ 2/3 r)}. 

In the following Proposition we summarize our main estimates formulated in terms of the 
scaled objects. The statement is a mere corollary of the previous lemmas. It follows by simple 
scaling from (|3TT^|) and (pT2U]) - (pr25)) . 

Proposition 1. There exists a constant C < oo, such that given any < r < r < oo and the 
scaled Lax entropy/flux pair defined as prescribed above, the following bounds hold uniformly in 



n : 



V'i(r_) D {(p,«):0<p<r-C- 1 ^H} 
T>%{r,r) C {(p,u) : p <r + CVr\u\} 



\s n PP { P , u ) 

\S n pu {p,u) 

\Suu(P, u ) 



< Cl.Tyn^{p,u) 



< 



< 



< 



< 



log (r/r) 
C 1 



log (r/r) r + p 3 
C 1 



C 



Ufa (r,r){P> u )> 
TLv%(r,r)(P, u), 

^V%(r,r)(P, u ), 



log (r/r) 3 



Rv$(r,r)(p,U) 



\F n {p,u)-^{p,u)S n Ap,u)\ < 



Cy/r 2 

-(r + u )lOT fcF) (p,it). 



(5.29) 
(5.30) 

(5.31) 

(5.32) 

(5.33) 

(5.34) 

(5.35) 

(5.36) 



P r ' n log(r/r) 
Due to (|5.29|) we can choose r large enough to ensure that for all n 

{(p,u) : pV\u\ < M} C P?(r), (5.37) 

where M is specified by the large deviation bounds given in Proposition |2] (via Lemma EJ. 
Further on, we choose r so large that 

-l 

(5.38) 



C 



< 100 sup \\og p(t,x)\ 
log(r/r) y (t, x )e[o,T]xT J 

and thus the bounds (|5.33|) - (|5.35j) of Proposition ^ are sufficient for our further purposes. 
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5.3 Outline of the further steps of proof 



In section El we present the main probabilistic technical ingredients of the forthcoming proof. 
These are variants of entropy inequalities and of the celebrated one and two block estimates. 
In Sectional we give an estimate for the terms with 'large' values of (p, u), we prove that 



Wo j GJn 



(5.39) 



< \ h n (t) + \ s n (t) + C [ h n (s) ds + o(l) 
* * Jo 



In Section |H1 we estimate the terms with 'small' values of (p,u), the section is divided into four 
subsections. 

In subsection 18. II we prove 



E 



1 

n 



I- 

n ' 



{ (M + pu n - pu) J n (p", u n )}(s, 

< Ch n (s) + o(l) 



(5.40) 



In subsection 18.21 we prove the one block estimate 



E 



1 

n 



{ (fa) n^(r - *(5f , c)) nr, ^)} («, 3 -) ds 



In subsection 18.31 we control the Taylor approximation 



(5.41) 



o 1 



E 



1 

n 



£ { (M (* n (r,u n ) - Tu n ) I n {T,u n )}(s, J -) 

< Ch n (s) + o(l) 



(5.42) 



Finally, in subsection 18.41 we control the fluctuations 



E 



~ E { G 5 " - p) ^ n - u ) r ^' (°> i) 



jGT" 



(5.43) 



< Ch n {s) + o{\). 



Having all these done, from and the bounds (fOHjl . (I5~1H|) . (I5~1T|) . (fOSJl . (I5~1H|) it 

follows that 



I A n (s)ds < - h n (t) + - S n {t) + C I h n {s)ds + o(l). 
io 2 2 y 



(5.44) 



Finally, from (|4.21j) . (|4.24j) . (|5.44j) and noting that s n (t) > we get the desired Gronwall 
inequality (|4.4|) and the Theorem follows. Note the importance of the term —dts n (t) on the 
right hand side of (|4.6|) . 
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6 Tools 

6.1 Fixed time estimates 

In the estimates with fixed time s £ [0, T] we shall use the notation 

L = L(n) := n~ 2f5 l. (6.1) 

Note that L 3> 1 as n — > oo. 

The following general entropy estimate will be exploited all over: 

Lemma 5. (Fixed time entropy inequality) 

Let I < n, V : fl l — > R and denote Vj(uS) := V(u>j, . . . , Wj+^i). Then for any 7 > 

E (^ E WW) < l -h-(s) + ±- ± Y, logE^ (exp {7Z, V,}) . 
jeT n ' ' jr'eT n 



(6.2) 



This lemma is standard tool in the context of relative entropy method. For its proof we refer 
the reader to the original paper |26| or the monograph |l(Jj . 

Proposition 2. (Fixed time large deviation bounds) 

(i) For any £ > there exists M < 00 such that for any s € [0, T] 

E G E {{ 1 +P n + \^\) 1 {^\u-\>M}}(s, J -)) <eh n (s) + o(l). (6.3) 

jeT n 

(ii) There exist C < 00 and M < 00 such that for any s 6 [0, T] 

E (^ E {l Sn | 2l {p"v|^|>M}}(^^)) <Ch n (s) + o(l). (6.4) 

The proof of Proposition's postponed to subsection llU.il It relies on the entropy inequality 
1)6. 2J) of Lemma El the stochastic dominations formulated in Lemma |H1 (see subsection 110. 1)1 and 
standard large deviation bounds. 

We shall refer to ()6.3)) and 1)6.4)1 as large deviation bounds. 

Proposition 3. (Fixed time fluctuation bounds) 

For any M < 00 there exists a C < 00 such that the following bounds hold: 



E(- 7" IS" - u\ 2 (s, < C h n (s) + o(l), (6.5) 
E (l E { I* 5 ™ " P\ 2 Mp-<m}}{s, J -)) < Ch n {s) + o(l). (6.6) 



The proof of Proposition|3]is postponed to subsection llU.21 It relies on the entropy inequality 
1)6.2)) of Lemma and Gaussian fluctuation estimates. 
We shall refer to ()6.5)) and 1)6.6)1 as fluctuation bounds. 
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6.2 Convergence to local equilibrium and a priori bounds 

The hydrodynamic limit relies on macroscopically fast convergence to (local) equilibrium in 
blocks of mesoscopic size Z. Fix the block size Z and (N, Z) E N x (wq/2)Z with the restriction 
N E [0, Z max 77], Z E [Z min £, Z max Q and denote 

I 1 

n l NjZ -.= {uen 1 ■. J2vj = N,J2Cj = z}, 

3=1 3=1 

I I 
K l NtZ (ui) := i^ fi {u | ^ Vj = N. ^ Cj = Z), 

3=1 3=1 

Expectation with respect to the measure n l N z is denoted by z ( ■ ) . For / : Q l N z — ► M. let 

1-1 

K l N ,zf(ui) :=J2zZ s(^3+i^',"%f(®t3 ! +i^ ~ /(^)» 

j=l w'flj" 

D l N , z (f) :=^EE^ ( £ a ( Wi , Wi+1 ;a/,^)(/(e^^-/(^) 2 ] . 

In plain words: ^Jy^ is the hyperplane of configurations w 6 1)' with fixed values of the 
conserved quantities, ir l N z is the microcanonical distribution on this hyperplane, K l N z is the 
symmetric infinitesimal generator restricted to the hyperplane Q l N z , and finally D l N z is the 
Dirichlet form associated to Kir z . Note, that K N z is defined with free boundary conditions. 

The convergence to local equilibrium is quantitatively controlled by the following uniform 
logarithmic Sobolev estimate, assumed to hold: 

(I) Logarithmic Sobolev inequality: There exists a finite constant H such that for any I E N, 
(N,Z) E N x (u!o/2)Z with the restriction N E [0,Zmax77], Z E [Zmin£,Zmax£], and any 
h : Q l N z — > R4. with Ejy^(Zi) = 1 the following bound holds: 

E l N>z (h log h) < N Z 2 ^ z (v^) . (6.7) 

Remark: The uniform logarithmic Sobolev inequality Q6.7JI is expected to hold for a very 
wide range of locally finite interacting particle systems, though we do not know about a fully 
general proof. In j^Tj the logarithmic Sobolev inequality is proved for symmetric K-exclusion 
processes. This implies that (|6.7(l holds for the two lane models defined in subsection [21 In [5] 
Yau's method of proving logarithmic Sobolev inequality is applied and the logarithmic Sobolev 
inequality is stated for random stirring models with arbitrary number of colors. In particular, 
(|6.7j) follows for the {— 1, 0, +l}-model defined in subsection [21 

The following large deviation bound goes back to Varadhan [221 • See also the monographs 
QUI and 01. 



36 



Lemma 6. (Time-averaged entropy inequality, local equilibrium) 

Let I < n, V : Q l — ► M+ and denote Vj(u) := V(u>j, ■ ■ ■ ,ujj + i-\). Then for any 7 > 



E (sEjfv,w)*j< 



(6.8) 



H Z 3 / 2 n 1+3/3+<5 1 



: 7 

Remarks: (1) Since 

_ 1+3/3+5 

^V- =o(1) ' 

in order to apply efficiently Lemma one has to chose 7 = 7 (ra) so that 

E^z(exp{ 7 V})=0(l), 

uniformly in the block size I = l{n) G N, and in iV G [0, Im&xrj] and Z G [Zmin£, I max(]. 
(2) Assuming only uniform bound of size ~ (HZ 2 ) 1 on the spectral gap of K l N z (rather than 
the stronger logarithmic Sobolev inequality (|6.7|) 1 and using Rayleigh-Schrodinger perturbation 
(see Appendix 3 of ^01) we would get 



E 



H^ 3 ||V||oo „ / m ^E^(V) max^Var^V) 

2n l +3/ 3+«5 S W+*II^Hoo^ ||v|| + ||v|| 2 



which would not be sufficient for our needs. 

(3) The proof of the bound (|6,8|) explicitly relies on the logarithmic Sobolev inequality l|6.7|) . It 
appears in |2H] and it is reproduced in several places, see e.g. [I], We do not repeat it here. 

The main probabilistic ingredients of our proof are summarized in Proposition 0] which is 
consequence of Lemma |H1 These are variants of the celebrated one block estimate, respectively, 
two blocks estimate of Varadhan and co-authors. 

Proposition 4. (Time-averaged block replacement and gradient bounds) 

Given a local variable £ : £l m — > M there exists a constant C such that the following bounds hold: 
(i) 

E (jT E(rT, C) } (s, x)\ 2 dx tbj < C (s n (t) + o(l)) . (6.9) 

(n) 

B (l I \ d ^ n ^ s ' x ^\ 2dxds ) <Cn 1 -^- 5 (s n (t) + o(l)). (6.10) 
(Hi) Further on, if £ : $7 — ► R (that is: it depends on a single spin) and £(o>) = whenever 
r)(uj) = then the following stronger version of the gradient bound holds: 
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The proof of Proposition 2] is postponed to subsection I1U.3I It relies on the large deviation 
bound 1(6.8)1 and some elementary probability estimates stated in Lemma ^] (see subsection 

MM- 

We shall refer to ()6.9)1 . respectively, 1)6.10(1 and (|6.11|) as the block replacement bounds, 
respectively, the gradient bounds. 

We shall apply (|6.9|) to £ = </> and £ = ip. From (|6.1U|) it follows that 

-t 



E 



J j \d x u n (s, x)\ 2 dxds^j < Cn 1 " 



(t)+o(l)), 



I j j \d x pT(s, x)\ 2 dxds) < Cn l+l3 -\s n {t) +o(l)). 



Using ()6.11() the last bound is improved to 
' -* f \d x (P{s,x)\ 2 



E 



Jf 



dxds) <Cn 1 - (3 - 5 (s n {t) + o(l)). 



The bound ()6.11|) will also be applied to £ = k (see ()2.11|) and ((2.12)1 ) to get 

/"* f \ n2P& f n{S : X) \ 2 dxds\ <Cn l -P- 5 {s n (t)+o(l)). 



o h 



(6.12) 
(6.13) 

(6.14) 
(6.15) 



7 Control of the large values of (p,u): proof of ( 15.391) 

7.1 Preparations 

In the present section we prove ((5.39)1 . First we replace i X^eT n ' ' ' ^y / T • • • dx. Note that 
given a smooth function F : T — > R 



1 

n 

Hence it follows that 

-< i 



n 



F(x) dx 



1/2 



(7.1) 



E / - E {(M {n 3/3 r)J n (r,u n )}(s,^-)ds 



(7.2) 



E 



J {(d x v) {n^r) J"(^,2")}( S ,x): 



where A™ 3 is again a simple numerical error term: 



d x ip n (s,x) 



dxds)+\lB \ I / ^^(sj^T dxds 

/O JT 



+ A /E ( / / \d x u n (s,x)\ z dxds 



Oin^r 1 ) = o(l). 



(7.3) 
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In the last step we use the most straightforward gradient bound (|4.12|) . (Using the gradient 
bound Q6.1UJI we could obtain the much better upper bound 

\A% i \ = 0(n(- 1 - s+7 W 2 ) = o(n 5e, r 1 ), 

but we do not need this sharper estimate at this stage.) 

So, we have to prove that the first term on the right hand side of (|7.2|) is negligible. Recall 
that J n = Sp. We start with the application of the martingale identity: 

E (7 {{vS n (f?\^)}(t,x) - {vS n (r,u n )}(0,x) - J* {{d t v)S n (r,u n )}(s,x)ds} dx 

E ([I V ^ X>> ( nl+/3Ln5n (^( x )' 2n ( x ))) ( X s) dxds 



(7.4) 

+ E ( I' I v(s,x) (n 1+/3+5 K n S n {p^(x),u n (x)^j (X?)dxds 



7.2 The left hand side of Ol 

From (|B3T|) . (|S321), and (|B37|) . we conclude that 

\S U {p,u)\ <C(p+\u\) l{p V |«|>M}- 

Hence, using the large deviation bound ()6.3|) it follows that, by choosing M sufficiently large 
we obtain 



S n {p^,u n ){s,-) 
n 



<eh n {s) + o{l). 



E (n ^ 

jeT n 

Hence, applying again (|7.1|) we get 

1 /"* 

|l.h.s. of (Ell < -h n {t) + C / h n {s)ds + o{\). (7.5) 

Remark: Note that this is the point where M and thus the lower edge of the cutoff is fixed. 
Also note the importance of the factor 1/2 in fromt of h n (t) on the right hand side. 

7.3 The right hand side of (J7.4|) : first computations 

First we compute how the infinitesimal generators n 1+ ^L n and n l+ P +& K n act on the function 
S n {p n {x),u n (x)): 

n 1+l3 L n S n (p n (x),v n (x))= (7.6) 
n 1+/3+s K n S n (p n (x),u n (x)) = (7.7) 
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where A™ 4 (x) and j4™ 5 (x) are the following numerical error terms: 

A^(x) = A™ 4 (uj_,x) := n l+f3 ^ ^ r(uj, u j+1 ; J, J) x 

S"^) + yK^) " a ( TO 7 " 1 ))(C / " 0)) 



-S;(r(x),u n (x))^a'(^-)(r,' - Vj ) 
-^(/T(x),S"(x))^a'(^)(C'-0)}. 



{^(W + ^(a(^) - a ( " X 7" 1 ))(^ " 



-S n {r{x),v n {x)) 



,2/3 



-5 p "(^(x),^(x))^a'(^)(r ? '-r ? ,) 



,/3 



-^(p"(x),^(x))^a'(^)(C'-C,)} 
+S£(p*(aO,tF(x)) x 

n l+3/3+5 , nx _j 1 nx _j -> 

— /2 \ a ( — j — " K i) + y a (— — 



+5:(/3"(r E ),S"(x)) 



E( , ,nx — js 1 „,nx — -I 
( a (— j )Uj+i-Xj) + ya ( 1 )Xj\- 



These error terms are easily estimated: using the fact that the second partial derivatives of S n 
are uniformly bounded and £ and r\ are bounded, by simple Taylor expansion after tedious but 
otherwise straightforward computations we find: 

sup sup|A? 4 (w,x)| < Cn 1+3(3 r 2 = o(l), (7.8) 

weSJ n rrGT 

sup sup \A^{ui, x) | < Cn 1+5 P +S r 3 = o(l). (7.9) 
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No probabilistic arguments are involved in these bounds. The global (averaged and integrated) 
error introduced by these terms will be of the same order. 

Next we do some further transformations on the main terms coming from the right hand 
sides of (j7.6|) and (|7.7|) . Performing integrations by part, introducing the macroscopic fluxes 
and using (|5.28|) we obtain: 

- J v(x)[s n p (T, u n )d x {n^r) + S:(fT, u n )d x (n 2 ^™) } (x) dx = 
J d x v{x) { {n^r) S^(p a , u n ) j (x) dx 

+ J d x v(x)[F n (r,u n ) - ^(F, u n )s n p {r, e^c*) dx 
+ J d x v{x)[ n ^s:(fP, u n ) (r - $(rr, a) } o«o <** 

+ / v(x){ n 3 ^ p (p n ,u n )(5 a! p n )(^-*(^,C h )) 

w i3 s n pu (r,u n ){d x u n ){r - * or\o) 
+n 2 X(/^' (^) - c n )) 

W^l(^,u n )(^)(^-$(5f,C n ))}(^)^ 

Note that, since J n = S*™, the first term on the right hand side is exactly the expression in the 
main term on the right hand side of 1)7.2(1 . Estimating the other terms on the right hand side 
of H7.1UJI is the object of the next subsection. 

Now we turn to the main term on the right hand side of 1)7.7(1 . Here, straightforward 
integration by parts yields 

-^v(x){s^r,un{n^diK n ) + s:(r,un{^dir)}(x)dx= (7.10) 

J d x v{x)[s n p (r, u n ) (n^d^) + SZ(r,u n ) {^d x x n ) } (x) dx 

+ I <x){ S; p (r,u n ){d x r){n^d x -k n ) 

+ S; u (r,u n ) ({d x u n ) {n 2 ?d x K n ) + {d x p n ) {n P d x T)) 

+ SUfP, u n ) {d x u n ) (/a,f) } (x) dx 
We will estimate the terms emerging from the right hand side in the next subsection. 
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7.4 The right hand side of (J7.4|) : bounds 
7.4.1 



We note that 

\F«(r, g") - ^(r, u n )s^r, u n )\ < c(i + ki 2 ) i { ^ v| ^ l>Af} , 

see ()5.36j) and 1)5.37(1 . Hence, applying the large deviation bounds (|6.3j) and Q6.4|) we obtain 



E 



da; 

< C/i n (s) + o(l). 



(7.11) 



7.4.2 

We use 



\s:(r,u n )\<c, 

see (|5.32|) and the first block replacement bound (|6.9|) to obtain: 



E 



JJ 



dx ds 



(7.12) 



<Cln^ l ~ 5+ W 2 = o{\). 



7.4.3 

Next we use 



sur,u n )\< 



c 



log (r/r) r + p 



\S n uu {fP ,u n )\ < 



sur,u n )\< 



c 



c 



log(r/r) + a/p" 



(7.13) 



log (r/r) ' 

see (|5.33[) . ()5.34|l . respectively, (|5.35jl . and note that here we do not exploit the fact that the 
constant factors on the right hand side are actually small. These, together with the block 
replacement bounds (|6.9|) . the gradient bounds (|6.12j) . (|6.14|) and the bound (|4.1() on the relative 
entropy s n (t) yield the following four estimates: 

{ n ^s; p (r,u^{d x r){r-^(r,c))}(s,x)\ dxds\<ci n ^ 5 = o(i), 



E 



E 



E 



E 



o JJ 
t 



JJ 

t 



JJ 

t 



{n^s; u (r,u n ){d x u n ){r-^(r,c n ))}(s,x)\ dxdsj <cin^- s = o(i), 
[n^s: P (r,u n )(d x r){d n -Hrx n ))}(s,x)\ d X d S ] <ci n - s = (i), p.u) 



10 JT 



{n^SZ u (r,^){d^){r-^,C))}(s,x)\ dxds) <Cln~ s = o(l) 
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7.4.4 



Using 



\s™(p n ,u n )\<c, \s:(r,u n )\<c, 



see (|5.31j) and (|5.32|) . and the gradient bounds (|6.12|) . I|6.13JI we obtain the following two bounds: 



n 



-1+/3+6 



E 



jt 



dx ds 

<Cn(- 1+w «/ 2 = o(l), 



n 



-1+/3+6 



E 



JT 



S. 



(7.15) 



dx ds 

<Cn(- 1 W)/2 =o(1) . 



7.4.5 



The following bounds are of paramount importance and they are sharp. We use ([7.13ft again 
and note that here we exploit it in its full power: the constant factor on the right hand side is 
small. These and the gradient bounds (|6.12j) and (|6.14j) yield the following three bounds: 



n 



n 



n 



n 



-1+/3+6 



E 



JT 



-1+/3+6 



E 



JT 



{s; p (r,u n ){d x r){n 2P d x 2 n )}(s,x) 
[s n pu {fP,u n ) (d x u n ) {n^d x K n ) } (a, x 



dx ds 
<cs n {t) + o(l), 



dx ds 
<cs n {t) + o(l), 



-1+/3+6 



E 



JT 



[s n pu {r,n n ){d x T){n P d x x n )\ 



s,x 



-l+(3+6 



E 



JT 



S^ u (r,u n ){d x u n ){^d x x n )}(s,. 



dx ds 
<c»"(t) 

ixis ) 

< cs n {t) + o{l). 



(7.16) 



The ratio r/r is chosen so large that 



c sup \v(t, x)\ < 
(t,x)e[0,T]xT 



1 



(7.17) 
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7.5 Sumup 

The identities lfTH|) . (|7Tjl . {TITty . (f7~TT!|l and the bounds JEBJ), f^SJ, dEDJ, f7?Tfl . (f7T4l . 
JUSJ), (TTTHl) yield 



E Uo J ^ dxV ^ ( n3 ^ n ) s p^i qn }}( s > x } dxds ) ~ ( r - h - s - of <El' 

1 /■* 

<-s n (t) + c h n (s)ds + o{l). (7.18) 

2 Jo 

Finally, from JOJ), JEHJ), (HI, JUJ) and (f7~T8|) we obtain 

8 Control of the small values of (p, u): proof of the bounds ( 15.401 ) 

to (GOES) 



8.1 Proof of rfPHj) 

We exploit the straightforward inequality 

\J n (r,U n )\ = 1^,2*) | < Cl {?nV | S n| >M}) 

see (|5.31[l and (j5.37j) . and boundedness of the functions p(t, x), u(t, x), d x v(t, x). Thus, applying 
the large deviation bound ()6.3|) we readily obtain (|5.4Uj) . 

8.2 Proof of (ICTTjl 

This is very similar to what has been done in various parts of subsection 17.41 We use the block 
replacement bound ()6.9|) and the bound 

\r{p n ,u n )\ = \l-S™(p? l ,u n )\<C (8.1) 

which follows from l|5.31j) . We readily obtain 



E 



"\^ n -V(v n ,C)) |/™(p n ,2 n )|}(s,x)| dsdx 

<C/n(- 1 - 5+3 «/ 2 = (l), 

which proves (|5,41|) , 

8.3 Proof of (I5TI21 

We write 

I n (fp,u n ) = l{p«v|u"|<Af} + ^{^v\u n \>M}I n (p n ,u n ), (8.2) 
and note that, by Taylor expansion of the function (p, u) i— > ^(p, u) 

\* n (F,vr)-Fv»\ l { ^v|2»|<M} < On" 2 ?. 
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On the other hand 



\m n (^,u n )\ < Cp^ 1 \u n \ 

and 

r \i n (r,u n )\ <c(i + |s n D, (8.3) 

see (E3S1) and (toTHI . Thus 

\V n (p n ,u n ) - p n u n \ \I n (p n ,u n )\ < C (n~ 2 ? + ( \u n \ + \u n \ 2 )l rp n Vl ^ l>M} ). 
From this, using the large deviation bounds (|6.3|) and (|6.4j) we obtain (|5.42|) . 
8.4 Proof of 

We use again (|8.2|) and (|8.3|) and get 

\{F-P) - «) I n (^, U n )\ < \{p n ~p) - U) I l { ^v|e»|<M} 

+C f 1 + In"! + |£ n | 2 'l 1 



{p n V|« n |>M} 

Now the fluctuation bounds (|6,5|) . (|6.6|) . and the large deviation bounds 1)6.3(1 . (|6.4|) together 
yield (JOE). 



9 Construction of the cutoff function: proofs 
9.1 Proof of Lemma El 

Proof. We sketch the proof for u > and leave the very similar u < case for the reader. Let 
pi > 0, u\ > be so chosen that for (p,u) £ [0, pi] x [0, u\] the following bounds hold with a 
fixed c > 0: 

|$ u (p, u) - * p (p, u) - (27 - l)u| < cu(u 2 + p), 
|$p(p,«) - 1| < c(n 2 +p), 

.2 



and 



|* u (p,M) -p\< cp(u + p), 



$ u (p, u) - * p (p, u)^0 for (p, «) / (0, 0). 



This can be done due to the (p,u) — ► (0,0) asymptotics of the macroscopic fluxes $ and Vf. 
It follows that as long as (a(u;r),u) £ [0, pi] x [0, iii] 

dp < 2p(l + c / (p + n 2 )) 

du ~ ^(2 7 - l) 2 u 2 + 4p + (2 7 - l)u" 

45 



This implies 

47-3 1 



/ 1- -t 7 -3 1 \ 
a(u;r) < r + C i \fru A r 4 t~ 2 u ^i- 1 I 



with a positive C, as long as (a(u;r),u) E [0,pi] x [0, u±\. 

Prom our assumptions it also follows that for p < pi and u > U\ 

dp , , 

where 

:= sup — -, — - < 00. 

p<pi P {$u{p,u) - ^ p {p,u)) 

U>U1 

Hence it follows that for u > u\ 

a(u; r) < a{u\;r) exp{6(n — ui)} 

as long as a(u;r) < p\. 

Putting these two arguments together the upper bound 

/ 1— 47-3 1 \ 
cr(u; r) < r + C\ I \fru A r 4 t- 2 u 2 ^- 1 I , u > 0, r < ro, 

follows with 

ro := sup |r : ^r + C ^/rui A r 4 T- 2 u^ 7 " 1 ^ ^ exp{6(n* — u\)} < p\\ , 

and 

r^(exp{b(u* -m)} - l) r , „ 

Ci = -2 V P1 V — '- + Cexp{6 n* - m)}. 

9.2 Proof of Lemma |H1 

Note first that given the bounds IJH.fJj) and condition ([Hj of subsection 13 . 1 1 1)5.24)) follows directly 
from (EH3I), and (|^TTT|) . So, we shall concentrate on ipT20]) - ()535)l only. 

By differentiating in the pde 1)5.1 1J) and applying straightforward transformations we obtain 
the following differential equations for S p , S u , S pp , and S pu , respectively: 

».(»«.)„> + (*„ - 9 P )(.S„)^ - * P (S P )„„ + (9.1) 

\ p / p \ p / p 

+ (*« - - * P (S U ) UU + (9.2) 



□ 
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*u(S pp ) pp + (*„ - * p ) (S pp)pu - 3> P (S pp ) uu + (9.3) 
+**>{lr) V„) p + 2*J?^*) (S pp ) u -Z p (f) S pp 

\ p J o \ P ' o \ P / oo 



p / p \ / p \^p / p p 

_2 ( ^lz!p \ s 

p \ % )pp pu ' 



*u(S pu ) pp + (*„ - ¥ p ) (Sp^ - *p(Spu) nu + (9.4) 



Because of the conditions outlined in subsection 13. II all the coefficients are smooth functions, if 
p is small enough. 

The respective initial conditions are 

S P (P^) = j^j^Wm + hf>e[r,r 0] }> (5 P )>,0) = 0, (9.5) 

S U ( P ,0)=0, (^)> > 0) = i ^^^ y l {pefc , r]} , (9.6) 

S PpM = ]^j^ 1 {pektf}> (S pp ) u ( P ,0)=0, (9.7) 
S pti (p,0) =0, 



Observe, that because of the asymptotics (|3.9|) we have 



1 + 0(P), -TT 1 ^ (p,0)=O(l). (9-9) 



p 



In order to understand the pdes (|9.1|) - ([9.4|) first we analyze in general the pde 

^uf PP + ($„ ~ ^p)f P u - $pfuu + Af p + Bf u + Gf = 0, (9.10) 

the functions {A, B, G} = {A, B , G}(p,u) being given on the left hand side of the pdes 
(|9.1|) - ([9.4|) . It is easy to check that A and G are even, B is odd with respect to u and also that 
k := A(0, 0) is 1, 27 — 1, 2, 27, respectively, in the four cases. 
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We solve in T> the Cauchy problem (|9.10f) with the initial condition 

f(p,0)=s(p), f u (p,0)=t( P ), P e[0,r ]. (9.11) 

The functions s(p) and t(p) will be identified with the various expressions in ()9.5|) - (|9.8j) . Then, 
in T>3(r) \ T> we solve the Goursat problem (|9.10j) with boundary conditions 

/(n,o) on dv 3 (f)\v, 

given by the solution of ~ _ ._. 

f y . on dVDV 3 (r). 

, the previous Caucny problem 



f(p,u) 



(9.12) 



Mind that /(n, 0) = 1 in the case of (jEEJ) and f(n, 0) = in the cases (j9~2l - (j9~4l . 
The pde (|9.10f) is hyperbolic in the domains considered. Its Jacobian matrix is 



D = J D(p,u):= . (9.13) 



The eigenvalues of D(p,u) are 



| = ±i | ^(<S> U - * p ) 2 + 4$ p * u T ($ u - * P ) } + $u (9.14) 



Mind that from the Onsager relation (|2.14j) it follows that for any (p,u) G "D ($ u — ^ p ) 2 + 
4$ p # u > 0. 

The characteristic coordinates (or Riemann invariants) w = w(p,u), z = z(p,u) of the pde 
(|9.10|) are determined, up to a functional relation 

w(p, u) = g{w(p, it)), z(p, u) = h(z(p, u)), (9.15) 

by the eigenvalue equations 

(w p ,w u )D = X(w p ,w u ), (z p ,z u )D = p(z p ,z u ). (9.16) 

Due to the gauge invariance (|9.15|) we can choose the characteristic coordinates w and z so that 

w(0, u) = ul {u>0} , z(0, u) = -ul {u<0} . 

This choice determines uniquely the characteristic coordinates. Observe, that as a corollary of 
Lemma |21 we have 

ci(Vp + u ) < w(p,u) < c 2 (VP + u) , (9-17) 

4 7 -3 1 

p < cz(p,u) 2 i- 1 «;(/), u) 2 ^- 1 . (9.18) 

w° := w(ro, 0) = z(tq, 0) =: z° 
w := w(r, 0) = z(r, 0) =: z 
w := w(r, 0) = z(r, 0) =: z 



We denote 
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In characteristic coordinates 

V = [0,u; ] x [0,z°], 

v n {u = 0} = [o, w°] x [o, z°] n{w = z}, 

V 3 {r)\V= [w°,u*] x [0,z] U [0,w] x [z°,u*], 
dV 3 (r)\V = {{w,z) : w G K,<u*]} U {(^,z) : z G [z°,u*]} 
dVr\V 3 (r) = {{w°,z) : z G [0, z°]} U {(w, z°) : to G [0,w°]}. 
The pde (j9,10j) written in characteristic coordinates reads 

f wz + af w + (3f z + uf = 0, (9.19) 

where 

X z - Au z + Bp z 



a = a(w, z) :-- 



A — a 



a of \ fJ'w Au w + Bp w 

(3 = P(w,z) := (9.20) 

A — p 

G(p w u z - p z u w ) 

v = u(w, z) := 

A — p, 

(Now all the functions on the right are understood as functions of (w,z).) In characteristic 
coordinates the initial conditions of the Cauchy problem in T> are 



(9.21; 



f(v,v) = s(p(v,v)) =: s(v), 

{fw ~ fz)(v,v) = 2u w {v,v)t(p(v,v)) =: t(v). 
The Cauchy problem ()9.19j) +( |9.21[) in the domain T> Pi {z < w} is solved by 

f(w ,z ) = -ip(zo,zo)s(zo) + -z(p{w ,wo)s(w ) (9.22) 

^ rwQ ^ y f w o 

+- / cp(v,v)t(v)dv + - / (tp w - ip z ) (v, v)s(v)dv 



9 ' ' 2 , 



■W{) 



+ / (p(v, v )(P(v , v) — a(v, v))s(v)dv 

J ZQ 

where the Riemann function (p is a solution of the adjoint Goursat problem: 

ip wz - (a(p) w - {f3<p) z + v<p = 0, (9.23) 

with boundary conditions: 

[ p(wo,t) = exp /* a(w ,v)dv, t G [z ,w ], 



Jz 

<p(s,z ) = exp f^ o (3(v,z )dv, sG [z ,w ]. 
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(9.24) 



Actually, the Riemann function depends also on (wo, zq): (p(w, z) = (p(wo, zo;w, z). In order to 
avoid heavy typography we omit explicit notation of this dependence. Note that in our cases 
(because of the left-right reflection symmetry of the respective pde), we will have (3{v,v) = 
a(v, v), thus on the right hand side of 1)9.22(1 the last term cancels. 
Also, if we consider the non-homogeneous pde 

h wz + ah w + (3h z + uh = g 

with the same initial conditions then it is solved by 

h(w Q ,z ) = f(w ,z ) + / / g(s,t)ip(s,t)dsdt, (9.25) 



where A WQtZ0 is the triangle with vertices (zq, zq), (wo,wo), (wo, zq). 

For details see any advanced textbook on partial differential equations, e.g. [S], 0, [3], HOI- 
In order to estimate / we will give a uniform estimate on the Riemann function (p. 

Proposition 5 (Bounds on the Riemann function). Let ip be the Riemann function asso- 
ciated to the equation 19.10\) (where the coefficients A,B,G are given on the left hand side of 
\9.1\) . i9,.^) or j9.4\ )) and wo > zq > 0. Then the following bounds hold uniformly for 

(wo, zq) S V with zo < t < s < wq and (s, t) £ V: 

K-l 

\<p(s,t)\ < c(— V 7 " 1 (9.26) 
I f s \^FT _1 



\(d w <p - d z <p)(s,s)\ < c 

w \W J 

Using Proposition with (|9.22|) . the initial conditions l|9.5|) - (|9.8j) and with Lemma we can 
estimate / in D which gives (|5.2(Jf) and (|5.21|) in this domain. 

The equations (|9.3|) and (|9.4j) are not closed for S pp and S pu , respectively, with the previous 
method one can only prove the required estimates for the solution of the respective homogeneous 
pdes. However, with (|9,25|) it is easy to show that these estimates can be extended for S pp and 
S pu , too. 

As an example, we show how to get the bound 1)5.22(1 for the homogeneous solution / of 
(|9.Hj) with initial conditions (|9.7|) . 

For the initial conditions we have the following bounds (using (|9. 17(1 . (|9. 18(1 ): 

c 

\s(v)\ < - — 2 ^{ve[z,z}}, 
log[r/r)v^ 1 l - 11 

\t(v)\ = 0. 
From Proposition El we have that for any zq < w$ 



i i 

V \ 2 T-1 ln / s n si C I V \ 2 T-1 



\(p(v,v)\<c[ — , \d w (p(v, v) - d z <p(v,v)\ < 

wo J w \w 



-1 
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Together with we get 



c 



i 



log (r/r) 

for any zq < wq. This means 



\f(w ,z )\ < ,L,, w 27 x i ^-^{zoe&zlzKwo}, (9-27) 



\f(w ,z )\ < ^ ^l^e^^wo} , 

which (using (|9.17|) ) translates to the (p, rz) coordinates the following way: 

c 1 

|/(p,u)| < : ^7^-lD 3 (r,r)(p,«)- 

log (r/r ) r K - ' 

(We only get this for u > 0, but from the symmetry of the pde this is also true for u < 0.) Also 
from (j9.27j) we get 

c L_ 2 L_ 



which (using (|9.18jl ) gives 



c 1 

|/0,«) I < ; 7=TT^-^V a (r.r){p,u). 

log (r/r ) p v - ; 



Putting together the two bounds on / we get the required estimate of (|5.22|) . 

Now back to the proof of Proposition |SJ The following lemma will be a basic tool for our 
estimates: 

Lemma 7 (Goursat estimate). Suppose the functions A(w , z) , B (w , z) , C (w , z) are defined 
on \x\,X2[ x [2/1,2/2] where < y% < 1/2 < x i < x 2 < 00 and /iaue f/ie following properties: 

sup / |-B(s,2;)|ds < -, 
zefoi>ie] •'ai " 

/•J/2 ^ 

sup / |A(u>,t)|dt < -, (9.28) 

\\xi ,Xl] J VI " 



r>IE 2 /-J/2 ^ 

/ |C(s,t)|<fedt < 

'i'l J 3/1 



6 



Let U be a solution of 



U wz - (AU) W - (BU) Z + CU = (9.29) 
on the rectangle [x\,X2] x [j/1,2/2] o,nd let 

sup \U(x\,z)\+ sup \U(w, y 2 )\ = M. 

J/l<Z<J/2 X1<W<X2 



Then 



sup \U(w,z)\ < 5M. 

(to,z)G[a:i,X2]x[j/i,j/2] 
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Proof. Denote U(w, 1/2) = f(w), U{x\,z) = g(z). Denote 

J(w) = f{w)- / B(s,y 2 )f(s)ds, 
g{z) = g{z) - f A(x u t)g(t)dt. 
Then U satisfies the following integral equation (for every (w,z) £ [£1,22] x [7/1,7/2]): 

U(w,z) = f(w) +g(z) -U(x 1 ,y 2 ) +/ A(w,t)U(w,t)dt + B(s, z)U(s, z)ds 



Jy 2 J xi 

C(s,t)U(s,t)dsdt. (9.30) 

'xi Jy-z 

Taking absolute values after some trivial estimates we get: 

IS 1 

\U(w,z)\<—M + - sup \U(s,t)\, 

D 1 {s,t)£[x 1 ,w]x[z,y 2 \ 

from that the needed bound follows immediately. □ 

Remark. Observe, that if we have some additional estimates on A Z ,B W and C then differenti- 
ating with respect to w or z and using the Gronwall inequality, one could also get bounds 
on \U W \ and \U Z \ in \x\,X2[ x [7/1,7/2]. 

Proof of Proposition 

(i) We first prove the proposition in the case when ^f(p, u) = pu and $>(p,u) = p + ju 2 . In that 
case the coefficient functions in (|9.1()|) take the following form: A = k, B = C = 0. Also, the 
Riemann invariants can be computed explicitly: 

W ( P ,u) := f VEE^t^A ^ ( V (27 _ 1W + 4p - ( 27 - 2)u) 



27-2 
47-3 



*l=i 27-2 

47-3 



(9.31) 



*M ■= ( V^g-^-D" ) ^ ( V (27 _ i) u 2 + Ap + (27 _ 2)M > 

One can easily check that the equations (|9.16j) hold. We define 

(3 {w) := (3(w,0) and U(w, z) := 99(77;, z) exp f- I (3 {s)ds) (9.32) 



((>() 



(for < z < w). Prom the definitions one can calculate that 

k 1 u w (w,0) « - 1 _! 
27 — 1 77,(77;, 0) 27 — 1 

Thus, it is enough to prove, that 

\U(s,t)\ < C, \d w U(s,s)\ < C-, |t9 2 *7(s,s)| < C (9.33) 

s 



52 



for zo < t < s < wo uniformly, with a constant depending only on k. To show (|9.33f) we will 
apply Lemma d 

From (SOU) and QKty we get that 



(aU) w -((P-(3 )U) z -a(3 U = 0, 



(9.34) 



and 



U(w ,z) 
U{w,z ) 



exp 



exp 



a(wo, t)dt 
(P(s,z )- 0(s,O))ds 



w 



Using the explicit formulas for the Riemann-invariants one can get estimates for the integrals 
needed for Lemma[7| Suppose [3:1,22] x [2/1,2/2] C [u,«;o] x [v,Zq] for some zo < v < wq. Then 
it can be shown that for zq < z < w < wq 



U-2 



a(w, t)dt 



rx2 



in 



< 



X 2 - X\ 



0(a, z) - 0(s,O)ds 



■•-1 



in 



(aflo)(s, t)ds dt 



< cz 



< 



w 

1 

x\ 



1 

X2 



(9.35) 



x\ x 2 . 



2/1 J 



where c only depends on k. From the first and second inequality it follows that \U(wo,z)\ and 
\U(w, zq)\ can be bounded by a constant depending only on k. Now fix s, t with zq < t < s < wo. 
Using ()9.35j) one can partition [v, wq] x [zq, v] into smaller rectangles in a way that the number of 
rectangles only depends on the value of k and on each small rectangle the conditions of Lemma 
13 hold (with A = a, B = (3 — (3q, C = —afto). Applying successively Lemma for the small 
rectangles (starting with the vertex (wo,zo)) one gets \U(s,t)\ < c(tz). Using the remark after 
Lemma one can also get the results of (|9.33|) for the partial derivatives of U. 
(ii) In the general case we do not know the explicit forms of the coefficients in (j9.1U|) only their 
asymptotics: 

A{p,u) = k(1 + 0{ P + u 2 )) 
B(p,u) = Cl u(l + 0(p + u 2 )) 
C(p,u) = c 2 {l + 0{p + u 2 )). 

We also do not have explicit formulas for the Riemann-invariants, but because ^(p, u) = pu(l + 
0(p + u 2 )) and $(p,u) = (p + ju 2 )(l + 0(p + u 2 )) if p, \u\ <C 1 the level-lines will approximate 
the respective level-lines of the system examined in (i). We will follow the steps of the proof 
for the specific case. We define (5q and U as in (|9.32j) . From the asymptotics we have 0q{w) = 
+0(1) which means it is again enough to prove ()9.33|) . We have the following equation 



-1 



27-I 

for U: 



(aU) w - ((/? - Pq)U) z + {u- a(3 Q )U = 0, 



(9.36) 
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with 



U(wq,z) = exp a(wo,t)dtj , 
U(w,z Q ) = exp I / ((3(s,z Q ) - (3(s,0))ds 



\Jw 

If we can prove similar bounds for the integrals of coefficients as in (|9.35|) then using Lemma 
the required estimates follow. From (|9.20|) we have that 

rx 2 ry2 r r Q 

I I u(s,t)dsdt = / / dpdu 

Jx! Jyi J J A A - I 1 

where A is the domain corresponding to the [xi,X2] x [1/1,1/2] rectangle on the (p,u) plane. If 
ro is small enough, then for (p, it) 6 T> 



C 



A — fl 



1 

< c 



\Ju 2 + p 



and since the right-hand side is integrable this gives uniform bounds on the previous integral. 
To get bounds on the integrals 

3/2 rx2 r%2 rV2 

a(w,t)dt, / (3(s, z) — /3(s, 0)ds, / / (af3o)(s,t)dsdt, 

Vi J xi Jx! Jyi 

we observe that because of the asymptotics described earlier if ro — > these integrals will be 
uniformly close to the respective integrals of the (i) case. Thus, again if ro is small enough (but 
fixed!) then the arguments of (i) may be repeated. □ 

We have proved that the bounds (|5.2U|) - (|5.23|I hold in T> if ro is small enough, now we have 
to extend this to the domain ^(r) \ T> for the solution of the respective Goursat-problems 
with boundary conditions (|9.12|) . The solution of these Goursat-problems can be expressed by 
integral equations similar to (|9.3fl|) (see |Hj , ) • Thus if we have some estimates for the solution 
on the boundary (which we have from (|9.12|) and the previous estimates of the Cauchy problem) 
then these may be extended (up to a constant multiplier) if we have uniform bounds on the 
integrals of the respective coefficients. 

If r is small enough then for (p, it) G D3 (r ) \ T> we have 

\$ P (p,u) - ® u (p,u)\ > c> 

which also implies that in this domain X — p > d > 0. From this it follows, that we can fix 
a small enough ro > such that in the domain T>s(fo) \ T> all the respective coefficients are 
well-defined smooth functions. Moreover, since this domain is compact, they are all bounded 
with a fixed constant which means that we have uniform bounds on the respective integrals. 
Using similar arguments as in the estimate of the solution of the Cauchy-problem one can get 
the bounds ()5.2Uj) - (|5.23|) also in this domain which completes the proof of Lemma |21 
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9.3 Proof of Lemma HI 

Proof. Note that 

(F-*S p ) p = * p S u -*S pp , 
(F-yS p ) u = $ u S u -VS pu . 
Also, there exists a constant C < oo such that for any (/?, u) G T> 

|*(p,«)| < cvH |$ P U«)| < c, |* u (p,u)| < c|«|. 

From these and the bounds (|5.21|1 . (|5.22)1 . (|5.23[) of Lemma |3] it follows that 

\(F - *S p ) p \(j>,u) < i -^( V ¥+ H) l V , 3[ r,r)( P ,u), 
\( F -* S p)v\M ^ log(^/r) (P+^\ U \)^ 3 (r_,r)(p,u). 

Integrating these and using (|5.16|) . the bound (|5.25|) follows. □ 

10 Proof of the "Tools" 

10.1 Proof of the large deviation bounds (Proposition 

Recall the definition (|fi.lj) of L. The following lemma follows from simple coupling arguments. 
Lemma 8. (Stochastic dominations) 

There exists a constant C depending only on max( s tX )e[0,T]xT p( s i x ) an d max (s,a;)e[o,T]xT \u(s, x)\ 
such that for any fixed (s, x) £ [0, T] x T the following stochastic dominations hold: 



P„n {p»(x) >z)<P (POI(L) > (z/C)L) , (10.1) 
P u J\u n (x)\ > z) < P(|GAU| > ((z/C) (10.2) 



where POI(L) is a Poissonian random variable with expectation L, and GAU is a standard 
Gaussian random variable. 

Lemma 9. (Large deviation bounds) 

(i) For any 7 < 00 there exists M < 00, such that for any n, j E T n and s £ [0, T] 

logE^(exp{ 7 L^)l { ^ )>M} }) < 1, 

logE„ ? (exp{ 7 L^)l {| ^ (i)|>M} }) < 1, 

(10.3) 

logE^(exp{7L|2"(^)|l {? „ (i)>M} }) < 1, 
logE„n(ex P { 7 L|^)|l {|2n( i )|>M} }) < 1. 
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(ii) For any 7 G (0, 1/(8C 2 )) there exists M < 00, such that for any n, j E T n and s £ [0, T] 

logE l/? (exp{7L|2"(^)| 2 l {?1(i)>M} }) < 1, 

(10.4) 

logE^(exp{7L|^(^)| 2 l {|SJ ^ )|>M} }) < 1. 

Proof, (i) We prove the first bound of (|1U.3|) . the other ones are done very similarly. 

Let Zl be a PO/(L)-distributed random variable. Using the stochastic domination (|10.1|) 
we obtain 



logE^(exp{ 7 L^)ll {?n( i )>M} }) 



<log(l + E w? (ex P { 7 L^)}l. 



< \/E, ? ( exp {2 7 Lj3»(£)}) J P, ? > M 



< JE{exp{iCZ L })Jp[(C/M)Z L > L 

< exp {|( (c (^_i) + (e 7C -!)-«)}, 

where a is arbitrary positive number. In the last step Markov's inequality is being used. Now, 
choosing a > exp(27(7) and M > (Ca)/(ln2) we obtain (|10.3j) . 

(ii) Again, we prove the first bound in Q1U.4|) . The other one is done in an identical way. 

Let again Zl be a PO/(L)-distributed and X be a standard Gaussian random variable. 
Using the stochastic dominations (jlU.lj) and (|1U.2|) we obtain 

logE ! , ? (exp{7L|2 n (^)| 2 l {?l(i)>M} }) 

<log(l + E^(exp{7L|n"(^)| 2 }l {? „ (i)>M} 



< \/E, ? ( exp {2 7 L|^)| 2 }) J P, ? (>(^) > M 



< yE(^exp{47C 2 (X 2 + L)}JyP^ L > (M/C)L 

< (1 - 8 7 C 2 )- 1/4 exp {^(4 7 C 2 + {e^' M - 1) - a)}, 

where a is arbitrary positive number. Given 7 < 1/(8C 2 ), we choose a sufficiently large and 
M > (Ca)/(ln2) to obtain fROjl . □ 

Now we turn to the proof of Proposition [2J 

Proof. The bounds ()6.3|) . respectively, ()6.4() follow directly from the entropy inequality (|6.2j) 
of Lemma El and the bounds (|l().3j) . respectively, (|1().4|) of Lemma 03 Recall that L ^> 1, as 
n — > 00. □ 
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10.2 Proof of the fluctuation bounds (Proposition OH) 

Within this proof we need the notation 

u n (s,x) : = ^^a(^^)(a-n-^( S ,^))=^(r E )-E^(u"(x)) : 



23 



,nx — k, 
1 I ' 



r/k-n 2/3 p(s, - 
v n' 



Since 



n 



<|i-tE«(^)II«(^)I + tE 



n 1 / 



u(s, ^) - u{s, -) 
n n 



I n 



o(l), 



and, similarly 



r(s,i.)- P (8±)\-\jr{8,±) 

\ fy-i / V ;V"J /I I V /v-| / 
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n 



n 



j — I I , j\ , 1 v /j — k s i / j \ / 



s|i-7E»(^) I !"(■.;) 



77 1 I 



n 



i T + - 

t 77 



o(l), 



we have to prove 



E (- V \u n (s,l)\ 2 ) <Ch n (s) + o(l), 



respectively, 



E (^ E I ^ 3 ~) I 2 h P~(s,i) i <m } ) * c h ^ s ) + 



(10.5) 



(10.6) 



Lemma 10. (i) There exists 7 > (sufficiently small) such that for all n, j £ T n and s £ [0, T] 

(10.7) 



logE,n(exp{ 7 L|5"( S ,^) | 2 }) < 1 



77' 



(ii) For any M < 00 there exists 7 > (sufficiently small) such that for all n, j 6 T n and 

se[0,T] 



log E,» ( exp { 7 L I \\ , ~ n(s i } , < M} } j < 1. 



(10.8) 
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Proof, (i) Let X be a standard Gaussian random variable, which is independent of all other 
random variables appearing in this paper, and denote by (. . . ) expectation with respect to X. 

logE^(exp{ 7 L|S"( S ,^) | 2 }) (10.9) 



= logE < (exp{^| ^(^-^(Cfc-E^CCfc)) | 2 

k 

= log (e^ ( exp {x^- J>(^) (a - E, ? (a: 

k 

Now, note that the random variables Oc — E l/ n(<^ fc ), k € T n , are uniformly bounded and under 
the distribution P u n they are independent and have zero mean. Hence there exists a finite 
constant C such that for any collection of real numbers Afc, k £ T n 

E^n ( exp { A fc (Cfc " E, ? (C fc )) }) < exp {C ^ Al} . 
Further on, there exists a finite constant C such that for any I 

y£Ky)| 2 <^ ( 10 - 10 ) 

k 

From these it follows that for some finite constant C, 

r.h.s. of (gnSD < log(exp{C7X 2 }> . 

Choosing 7 sufficiently small in this last inequality we obtain ()10.7j) . 

(ii) Note first that, given M < 00 fixed, there exists a zero mean bounded random variable 
Y such that for any r£K 



r 2 



l{|r|<A/} < logE( expjry 



Let Y\,Y2, . . . be i.i.d. copies of Y which are also independent of all other random variables 
appearing in this paper, and denote by (...) expectation with respect to these. Then we have 

log E u n ( exp { 7 L I r{s, 3 ~) \\ , ^ , < M} }) (10.11) 

< log (E^n ( exp { p= ^ P a (^—j— ) {Vk ~ E K (%) 

k 

Next note that for any A < 00 there exists a constant C < 00 such that for any n € N, any 
s G [0, T] and any collection of real numbers A& £ [— A, A], A; G T n 

E^ ( exp { ^ A fc (% - B u n (%))})< exp {Cn" 2 ^ ^ A 2 } . 
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Hence, using again 1)10.10)1 . 

r.h.s. of (HnHD <log(e X p{C7 + . . .+ Y^) / }). 

Now, since the i.i.d. random variables Yi,Y2,... are bounded and have zero mean, choosing 
7 sufficiently small this last expression can be made arbitrarily small, uniformly in L. Hence 

cnn). □ 

Now back to the proof of Proposition 

Proof. From 1)6. 2J) and 1)10.7)) . respectively, from 1)6.2)1 and 1)10.8)) we deduce ()10.5)) . respectively, 
()10.6)) . Finally, these two bounds and the arguments at the beginning of the present subsection 
imply 1)6.5)1 . respectively, 1)6.6)1 . □ 

10.3 Proof of the block replacement and gradient bounds (Proposition HJ) 
10.3.1 An elementary probability lemma 

Let (0, 7r) be a finite probability space and Wj, i 6 Z i.i.d. Si-valued random variables with 
distribution ir. Further on let 

£:ft m ^R, &:=£(w i ...,w i+ro _ 1 ). 
For x 6 co(Ran(£)) denote 

v _ E 7r (eiexp{E^ 1 A-U}) 
E^explA-Ci}})™ ' 

where co(-) stands for 'convex hull' and A £ M. d is chosen so that 

E^dexpjA-Ci}}) 
E^(exp{A-Ci}}) 

For leNwe denote plain block averages by 

Finally, let 6 : [0, 1] — > R be a fixed smooth function and denote 

M(b):= f b(s)ds, V(b) := M(b 2 ) — M(bf 
Jo 

We also define the block averages weighted by b 

l i 



x. 



< fe ' 0/ == yEHi/OC;, (6, 0i : = yEW. 

3=0 3=0 

The following lemma relies on elementary probability arguments: 
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Lemma 11. (Microcanonical exponential moments of block averages) 

There exists a constant C < oo, depending only on m, on the joint distribution of (£i,Cj) an d 
on the function b, such that the following bounds hold uniformly in I G N and x £ (Ran(£) + 
••• + Ran(C))A: 

(i) IfMib) = 0, then 

E(exp{ 7 v^<6,£)i}|C/ = *) <exp{C(7 2 + 7 M)}. (10.12) 

(ii) IfM(b) = 1 then 

E( exp { 7 v^((6 , Oi ~ 3((6 , C)i)) } \ Ci = x) < exp{C( 7 2 + y/Vt)}. (10.13) 

Proof. We prove the lemma with m = 1, that is with (£i)' =1 independent rather than in- 
dependent. The m-dependent case follows by applying Jensen's inequality in a rather straight- 
forward way. 

(i) In order to simplify the argument we make the assumption that the function s \— > b(s) is 
odd: 



6(1 -s) = -b(s). 



(10.14) 



The same argument works if the function s \— > b(s) can be rearranged (by permutation of finitely 
many subintervals of [0, 1]) into a piecewise continuous odd function. This case is sufficient for 
our purposes. The proof of the fully general case — which goes through induction on I — is 
more tedious and it is left as a fun exercise for the reader. 
Assuming (|10.14[) we have 

m\ 

v^<&,Oj = z- 1/2 £&0'/0&-&-i) 

3=0 

and hence 

E(exp{ 7 v / Z(6,0i}|0 = a ; ) 

= E(E(exp{ 7 v / /(6, Oi} Cj+Ci-j- J = 0, . . . , Z 



Cj + C 



1-3 



Ci = x 



[1/2] 

= e( h E(ex P {yi- i/2 bu/im - e»-i)} 

3=0 

[1/2] 

^expjc^^r^ovo 2 } 

3=1 

= exp {C 1 \V(b) + 0(1/0)}- 

In the second step we use the fact that the pairs (£j>£t-j)) j = 0, . . . , [1/2] are independent, 
given £ • + Ci-j> J = 0, . . . , [Z/2]. In the third step we note that the variables £j are bounded 
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andEfe-^ICi + Ci-i) =0. 

(ii) Beside S(a;) we also introduce the functions 

Ei : (Ran(C) + • • • + Ran(C))// R, 3, (as) := E(a|Cz = x). 
We shall exploit the following facts 

(1) The functions 3(sd) and 3^(a:) are uniformly bounded. This follows from the boundedness 

(2) The function x 3 (a;) is smooth with bounded first two derivatives. This follows from 
direct computations. 

(3) There exists a finite constant C, such that 

|3,(x)-H(a!)| <Cr\ 

This follows from the so-called equivalence of ensembles (see e.g. Appendix 2 of |1U|)- 
We write 

(b,0i-Z((b,Qi) = «6, Oi - Ci) + (?« - SiCCi)) (10-15) 

+ (SKCi)-3(C J )) + (3(C0-3((6,CW)- 

By applying Jensen's inequality we conclude that we have to bound the exponential moments 
of type (jlO.l.Sj) . separately for the four terms. 

Bounding the first and last terms reduces directly to (|10. 12|) . the third term is uniformly 
0(Z _1 ), so we only have to bound the exponential moments of the second term in (|10.15j) . This 
is done by induction on I. Let C(l) be the best constant such that for any 76M 

E(exp { 7 VZ(& -S,(? { ))} Ci = x) < exp{C(07 2 }- 



We prove that C(l) stays bounded as I — > 00. 
The following identity holds 



1 



v^TT 



(6+1 - 3i(C;+i)) 



^-py (3;(C;) - 3 ;+ i(Ci+i)) 

' s i(Cz+i) - s H-i(Ch-i)) 



VTTT 
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Thus 



E(exp{ 7 VlTT(e l+ i-S !+ i(C w ))} Ci+i =x 
E(E(exp { 7V TTT(e m - 3 I+1 (C,+i)) } 

7^ 



E 



E 



E 



(expj-^^-SKCO)} 



C« j Cm J Cz+i — * 



c 



i+1 



7/ 



exp{-^===(H z (Ci) - ~;+i(C m ))} 



exp{ 



7 



\/T+T 



l(C 1+1 ) 



(Ci+i))} 



X 



The terms 



fo+i - Si(Ch-i)), «( 3 »(Ci) - Sj+iCCi+O), (Si(Cj+i) - 3i +1 (Ci+i)) 
are uniformly bounded and 

E(?h.i-3i(C i+1 )|Ci+i) = 0, 
E(3,(C,)-S,+i(Ci +1 )|Ci+i) = o, 
E(S 1 (C J+1 )-Sj +1 (Cz +1 )|Cz+i) = o. 
Using the induction hypothesis it follows that there exists a finite constant B such that 

c(i+1) <J_c. (i) + _L_ B . 

Hence, lim sup^^ C(Z) < .B and the lemma follows. 



□ 



Lemma 12. (Microcanonical Gaussian bounds) 

There exists a 70 > 0, depending only on m, on the joint distribution of (£j,CJ an d on the 
function b, such that the following bounds hold uniformly in I £ N and x £ (Ran(£) + • • • + 
Ran(C))/Z: 

(i) IfM(b) = 0, tfien 

logE(exp{ 7 oZ(&,0/ 2 }|0 = z) <!• (10.16) 

(ii) IfM(b) = 1 then 



logE(exp{ 70 /((6,0z-3((6,C) i )) 2 } 



C, = a; < 1. 



(10.17) 

Proof. This is actually a corollary of Lemma 1111 The bounds 1)10.16(1 and (|1U.17() follow from 
(|10.12ft . respectively, <|10. by exponential Gaussian averaging (as in the proof of Lemma 

hdj. □ 
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10.3.2 Proof of Proposition SI 

Now we turn to the proof of Proposition 

Proof, (i) In order to prove (|6.9|) first note that by simple numerical approximation (no proba- 
bility bounds involved) 



c _ 



We apply Lemma El with 



V = \{e- E(rT, C)} (0) | 2 = | (a , 0i ~ H«a , r,) lt (a , t ) 



We use the bound (|10.17j) of Lemma H2l with the function b = a. Note that 7 = 70/ can be 
chosen in (|6.8|) . This yields the bound (|6,9I) , 

(ii) In order to prove (|6.10|) we start again with numerical approximation: 



/ \d x C(x)\ 2 dx-- V \d x C( 3 - 

We apply Lemma El with 



J N |2 



We use now the bound (|1U.16|) of Lemma ^] with the function b = a'. Now we can choose 
7 = 7ol 3 /n 2 and this will yield the bound (|6.10|) . 
(hi) Next we prove (|6.11j) . We apply Lemma El with 



V 



d x e(0)\ 2 = n 2 I ZkAk/l)Zk 
rjn(0) P Ek<k/l) Vk 



n 



2 I £ fc a'(V0(&-£-fc) 



2/3 Z k a ( k MVk+V-k) ' 

where in the last equality we use the fact that the weighting function x *— > a{x) is even. We 
compute the exponential moment E 2 ^ 1 ( exp{7V}) . Let X be a standard Gaussian random 
variable, which is independent of all other random variables appearing in this paper and denote 
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by (...) averaging with respect to it. We have 
EK(exp{7V} 



E 2l+1 ( exp It— I , 



E fc o(V0 (»/* + »/-* 

E^fexpjXVr 



Z 3 / 2 v / L#A)fe+^) 



< ^ { exp { 7 ^ ^ fjfc/n( — , } J ) 



<(exp{CX 2 7 ^}), 



i 3 

where we used the facts that the random variables % are non-negative, O is finite and rj(uj) = 
implies £(ui) = 0. In the last step we used the inequality 

a'(xf < Ca(x), 

which follows from the conditions on a(x), see subsection 14.31 

From this bound it follows that in Lemma we can choose 7 = 7oZ 3 /n 2 , with a small but 
fixed 70, and hence the second bound in (|6.11|) follows. 

□ 

11 Appendix: Some details about the PDE ( 11.11) 

Hyperbolicity: One has to analyze Jacobian the matrix 

D . = ( ( P ^P ( PU ^ \ = ( U 9 \ 

\ {p + ju 2 ) p {p + ju 2 ) u J \1 2 7 J 
The eigenvalues with the corresponding right and left eigenvectors are: 

Dr = Xr, Ds = fj,s, l^D = XI , nvD = /xm , 
(ft stands for the transpose of the column 2-vector v). The eigenvalues and eigenvectors are 
X 



/' 

and 



it 

t 



m 



\=±l { V(2J - l) 2 ^ 2 + 4/o ± (2 7 + 1)«, } 

l - { Tv / (27-l) 2 " 2 + 4/9 - (2 7 - l)u} , l) , 
1,~ {±^(27 " 1) 2 ^ 2 + 4p - (2 7 - 1)^}) • 
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So, we can conclude that the pde (jl.lj) is (strictly) hyperbolic in the domain 



7^1/2: 
7 = 1/2: 



{(p,u)eR+xR: (p,u) ^(0,0)} 
{(p,u) £l + xl:p/0}. 



Riemann invariants: The Riemann invariants w = w(p,u), z = z(p,u) of the pde are given by 
the relations 

(w p ,w u ) ■ s = = (z p ,z u ) ■ r. 

That is, the level lines w = const., respectively z = const, are determined by the ordinary 
differential equations 



Actually only the level lines of the functions w(p,u), respectively, z(p,u) are determined. In 
our case the Riemann invariants can be found explicitly. For 7 / 3/4 we get 



Where F : M. — > M. is an appropriately chosen bijection (mind, that only the level sets of the 
Riemann invariants are determined). 

Note that due to the changes of sign of 27 — 1 and 27 — 2, the above expression gives rise to 
qualitatively different behavior of the Riemann invariants. The picture changes qualitatively at 
the critical values 7 = 1/2, 7 = 3/4 and 7 = 1. In Figure |3 we present the qualitative picture 
of the level lines of w(p,u) and z(p,u) for 3/4 < 7 < 1, and 7 > 1, respectively. (For economy 
reasons we omit the graphical representation of the other cases, but encourage the reader to 
sketch it.) In all cases the Riemann invariants satisfy the convexity conditions 



in x M. for all 7. (The sign of the function F(-) is so chosen, that these expressions be non- 
negative.) The inequalities are strict in the interior of M+ x M, except for the 7 = 1 case, when 
these expressions identically vanish. These conditions are equivalent to saying that the level 
sets {(p, u) G [0, 00) x (— 00, 00) : w(p, u) < c} and {(p, u) £ [0, 00) x (—00, 00) : z(p, u) < c} be 
convex. See |llj . |12j or for the importance of these convexity conditions. 

It is of crucial importance for our problem that the level curves w(p,u) = w = const, 
expressed as mm p(u, w) are convex for 7 < 1, linear for 7 = 1 and concave for 7 > 1. 

Genuine nonlinearity: Genuine nonlinearity holds if and only if 




w(p,u) 
z(p,u) 





(11.1) 



(A p , X u ) • r / / (p P ,Pu) • s. 
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(a) 3/4 < 7 < 1 



(b) 1 < 7 



Figure 3: Level lines of Riemann- invariants 



in the interior of the domain R + x R. Elementary computations show that 



(A p , A u ) • r = 
(MpjMu) ■ s = 



P = - 



47(27 ~ I) 2 
(7 + I) 2 



u 



2 and 



{ 



u < 
u > 



(11.2) 



Thus, for 7 > 0, 7 / 0,1/2 the system is genuinely nonlinear on the closed domain R + x R; 
for 7 = 0, 1/2 it is genuinely nonlinear in the interior of R+ x R (with genuine nonlinearity 
marginally lost on the boundary, p = 0). For 7 < genuine nonlinearity is lost in the interior 



Lax entropies and entropy solutions: Lax entropies of the pde are solutions of the linear 
hyperbolic partial differential equation 

pS pp + (27 - l)uS pu - S uu = 0. 

It turns out that the system is sufficiently rich in Lax entropies. In particular a Lax entropy 
globally convex in R + x R is 



Construction of other Lax entropies with particular features (e.g. possessing scale similarity, or 



The Maximum Principle and positively invariant domains: For 7 > our systems satisfy the 
conditions of the Lax's Maximum Principle proved in jllj . Namely: (i) they do posses a globally 
strictly convex Lax entropy bounded from below, see (jll..3j) : (ii) the Riemann invariants w(p, u) 
and z(p, u) satisfy the convexity condition (iii) they are genuinely nonlinear in the interior 

of V, see (HSJ). 

Hence it follows that convex domains bounded by level curves of w(p,u) and z(p,u) are 
positively invariant for entropy solutions. 



of R , x R. 




(11.3) 



polynomial in ^fp and u, etc.) is a very instructive exercise. 
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First we conclude, that T> itself is positively invariant domain, as it should be. 

Second: a very essential difference between the cases 7 < 1, 7 = 1 and 7 > 1 follows, which 
is of crucial importance for the main result of the present paper. In the case 7 < 1 all convex 
domains bounded by level curves of the Riemann invariants are unbounded (non- compact) and 
thus there is no a priori bound on the solutions. Even starting with smooth initial data with 
compact support nothing prevents the entropy solutions to blow up indefinitely after appearence 
of the shocks. On the other hand, if 7 > 1 any compact subset of T> is contained in a compact 
convex domain bounded by level sets of the Riemann invariants, which fact yields a priori 
bounds on the entropy solutions, given bounded initial data. A microscopic consequence of this 
fact is that the proof of our main theorem is valid only for 7 > 0. 
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